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ABSTRACT 

Weakly collisional magnetized cosmic plasmas have a dynamical tendency to develop 
pressure anisotropies with respect to the local direction of the magnetic field. These 
anisotropics trigger plasma instabilities at scales just above the ion Larmor radius pi 
and much below the mean free path Amfp. They have growth rates of a fraction of 
the ion cyclotron frequency, which is much faster than cither the global dynamics or 
even local turbulence. Despite their microscopic nature, these instabilities dramati- 
cally modify the transport properties and, therefore, the macroscopic dynamics of the 
plasma. The nonlinear evolution of these instabilities is expected to drive pressure 
anisotropics towards marginal stability values, controlled by the plasma beta Pi. Here 
this nonlinear evolution is worked out in an ab initio kinetic calculation for the sim- 
plest analytically tractable example — the parallel {k± = 0) firehose instability in a 
high-beta plasma. An asymptotic theory is constructed, based on a particular physi- 
cal ordering and leading to a closed nonlinear equation for the firehose turbulence. In 
the nonlinear regime, both analytical theory and the numerical solution predict secu- 
lar (cx t) growth of magnetic fluctuations. The fluctuations develop a spectrum, 
extending from scales somewhat larger than pi to the maximum scale that grows sec- 
ularly with time (oc t^^^); the relative pressure anisotropy (jp±^ — tends to the 
marginal value — 2//?^. The marginal state is achieved via changes in the the magnetic 
field, not particle scattering. When a parallel ion heat flux is present, the parallel fire- 
hose mutates into the new gyrothermal instability (GTI), which continues to exist up to 
firehose-stable values of pressure anisotropy, which can be positive and are limited by 
the magnitude of the ion heat fiux. The nonlinear evolution of the GTI also features 
secular growth of magnetic fiuctuations, but the fiuctuation spectrum is eventually 
dominated by modes around a maximal scale ~ yoJr/Amfp, where It is the scale of 
the parallel temperature variation. Implications for momentum and heat transport are 
speculated about. This study is motivated by our interest in the dynamics of galaxy 
cluster plasmas (which are used as the main astrophysical example) , but its relevance 
to solar wind and accretion flow plasmas is also briefly discussed. 

Key words: galaxies: clusters: intracluster medium — instabilities — magnetic fields — 
MHD — plasmas — turbulence . 
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It has recently been realized in various astrophysics and 
space physics contexts that pressure anisotropies (with re- 
spect to the direction of the magnetic field) occur natu- 
rally and ubiquitously in magnetized weakly collisional plas- 
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mas0 They lead to very fast microscale instabilities, fire- 
hose, mirror, and others, whose presence is likely to funda- 
mentally affect the transport properties and, therefore, both 
small- and large-scale dynamics of astrophysical plasmas — 
most interesti ngly, the plasmas of galaxy clusters and ac- 
cretion discs ([Hall fc Sciamalll979l: ISchekochihin fc Cowle" 



20061 : ISchekochihin et all l2005f . l2008l: ISharma et alj |200l. 
2007l : lLvutikovl2007^ . These instabilities occur even (and es- 



pecially) in high-beta plasmas and even when the magnetic 
field is dynamically weak. The current state of theoretical 
understanding of this problem is such that we do not even 
have a set of well-posed macroscopic equations that govern 
the dynamics of a plasma in which the coUisional mean free 
path exceeds the ion Larmor radius, Amfp ^ pi (equivalently, 
ion collision frequency is smaller than the ion cyclotron fre- 
quency, Vii <^Q.i). This is because calculating the dynamics 
at long spatial scales I ^ pi and slow time scales correspond- 
ing to frequencies uj <Si 0,% requires knowledge of the form 
of the pressure tensor and the heat fluxes, which depend on 
the nonlinear evolution and saturation of the instabilities 
triggered by the pressure anisotropies and temperature gra- 
dients. Since this is not currently understood, we do not have 
an effective mean-field theory for the large-scale dynamics. 

In the absence of a microphysical theory, it is prob- 
ably sensible to assume that the instabilities will re- 
turn the pressure anisotropies to the marginal level 
and to model large-scale dy namics on th i s ba si s, via 



a suitable closure scheme (ISharma et al. |2006|. 20071: 
Schekochihin fc CowlevI l2006l : iLvutikovl 120071 : iKunz et al.l 



20111 ). Thi s approach appears to be supported by the solar 
wind data (IGarv et al."200l':'Kaspcr. Lazarus fc Gar^l2002l: 



Marsch. Ao fc Tul2004 : HcUingcr ct al,,2006 : .Matteini et all 



20071 : ISale et al ll2009l ). However, a first-principles calcula- 



tion of the nonlinear evolution of the instabilities remains 
a theoretical imperative because, in order to construct the 
correct closure, we must understand the mechanism whereby 
the instabilities control the pressure anisotropy: do they 
scatter particles? do they modify the structure of the mag- 
netic field? The calculation presented below will lead us to 
conclude that the latter mechanism is at work, at least in 
the simple case we are considering (see discussion in section 
16. If) , and indeed a sea of microscale magnetic fluctuations 
excited by the plasma instabilities will act to pin the plasma 
to marginal stability. 

In this paper, we present a theory of the nonlinear evolu- 
tion of the simplest of the pressure-anisotropy-driven insta- 
bilities, the parallel (k ± — 0) firehose instability and the gy- 
rothermal instability (|Schekochihin et al.ll2O10l ). To be spe- 
cific, we consider as our main application a plasma under 
physical conditions characteristic of galaxy clusters: weakly 
coUisional, fully ionized, magnetized and approximately (lo- 
cally) homogeneous. We will explain at the end the extent 
to which our results are likely to be useful in other contexts, 
e.g., accretion flows and the solar wind (section [T]). 

The plan of exposition is as follows. In section[2l we give 



^ As will be explained in detail in what follows, by weak coUision- 
ality we mean a state where Larmor motion is much faster than 
the collision rate, but large-scale dynamics occur on time scales 
slower than collisions, so collisions neither can be neglec ted nor 
are th ey sufhciontly dominant to justify a fluid closure. iBalbuj 
1I2OO4I ') calls this state a "dilute" plasma. 



an extended, qualitative, mostly low-analytical-intensity in- 
troduction to the problem, explain the relevant properties 
of the intracluster plasma (section 12. If) , the origin of the 
pressure anisotropies (section 12. 2p , sketch the linear theory 
of the firehose instability (section I2.3[l . the main principle 
of its nonlinear evolution (section 12. 4p . and show that a 
more complicated theory is necessary to work out the spa- 
tial structure of the resulting "firehose turbulence" (section 
I2.5|l . In section [31 a systematic such theory is developed via 
asymptotic expansions of the electron and ion kinetics (the 
basic structure of the theory is outlined in the main part 
of the paper, while the detailed derivation is relegated to 
Appendix|X]), culminating in a very simple one-dimensional 
equation for the nonlinear evolution of the firehose fluctu- 
ations (section 14. ip . the study of which is undertaken in 
section [l] The results are a theoretical prediction for the 
nonlinear evolution and spectrum of the firehose turbulence 
(section 14.31) and some tentative conclusions about its effect 
on the momentum transport (section 14. 4[) . In section [5l we 
extend this study to include the effect of parallel ion heat 
flux on the firehose turbulence: in the presence of a parallel 
ion temperature gradient, a new instabil ity emerges (the gy- 
rothermal instablity recently reported bv lSchekochihin et al.l 
I2OIOI and recapitulated in section — which, under some 
conditions, can take over from the firehose. For it as well, 
we develop a one-dimensional nonlinear equation (section 
I5.1[l . solve it to predict the nonlinear evolution and spatial 
structure of the gyrothermal turbulence (section 15.31) and 
discuss the implications for momentum transport (section 
15. 4|) . A discussion of our results and of the ways in which 
they differ from previous work on firehose instability in col- 
lisionless plasmas is given in section (6] A brief survey of 
astrophysical implications (both galaxy clusters and other 
contexts) follows in section [T] Finally, section [8] contains a 
very concise summary of our findings and of the outlook for 
future work. Note that while section [2] is largely a pedagogi- 
cal review of our earlier work ([Schekochihin fc Cowlevll2006l : 
ISchekochihin et al.ll2005l . |2008| ). most of the theorv and re- 
sults presented in sections [3HS1 is new. 

A reader not interested in the technicalities of kinetic 
theory is advised to ignore section [3] and Appendix [X] A 
reader only interested in the formal derivation may skip sec- 
tion [21 as section [31 (supplemented by Appendix IX)) and the 
sections that follow it can be read in a self-contained way. 



2 QUALITATIVE CONSIDERATIONS 

2.1 Galaxy clusters: observations, questions, 
parameters 

Galaxy clusters have long attracted the interest of both 
theoreticians and observers both as dynamical systems 
in their own right and as co smological probes (jBahcalll 
I2OOOI : [Peterson fc Fabianll2006l ). While gravitationally they 
are dominated by dark matter, most of their luminous 
matter is a hot, diffuse , fully ionized, X-ray emitting 
hydrogen plasma ([Sarazi nl l2003l ) known as the intracluster 
medium, or ICM (the galaxies themselves are negligible 
both in terms of their mass and the volume they occupy). 
Crudely, we can think of an observable galaxy cluster as 
an amorphous blob of ICM about 1 Mpc across, sitting 
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in a gravitational well, with a density profile peaking at 
the center and decaying outwards. Observationally, on 
the crudest level, we know what the over all density and 
temp er ature profiles i n clust er s are (e.g., Vikhlinin et al.l 



20051: ■ , , 

Cavagnolo et al.l l2009l l . Recent 
observations reveal the ICM to 
multiscale structure displaying 



Piffaretti et"all l2005l : iLeccardi fc Molendil l2008l : 



ments, waves 



highly resolved X-ray 
be a rich, complicated, 
r ipples, bubbles, fi la- 
shocks, edges etc. ([Fabian et al.l l2003a|| b| . 



2005al . [200^: ISanders fc FabianI [iood . l2008l : iForman et all 



2OO7I : iMarkevitch fc VikhMnl l2007h. temperature f luctu- 



ations JSimionescu et aL 200ljr" Markevitch et all l2003l: 



Fabian et all I2OO6I: iMillion fc AllenI l2009l: ISanders et all 
2010al : Lagana. Andrade-Santos fc Lima Netd 20101 ) 
and most proba bly also broad - band disordered tur- 
bulen t motions jChur azov et al.l '2004'; ' S chuccker et al 



2004 



2ooe 



Rebusco et al.l 



Sanders et al.l 



2005 
l2010bl . 



2006 
I2OIII : 

that 



2008; Grahara et al 
lOgrean et all |2010D 



Radio observations tell us 
tangled magne tic fields, which 



the ICM also hosts 
are p robably dynami- 



20041: Vogt fc Enfilin 20051: Kuchar fc Enfilinl 20091: 


Clarke fc Enfilin 2006: 


Govoni et all 20061: ICuidetti et all 


20081: iFerrari et al.ll2008 


). 



These and other observations motivate a number of 
questions about the ICM, which are representative of the 
problems generally posed for astrophysical plasma systems 

• Can we explain the observed ICM temperature pro- 
files, in particular the apparent lack of a cooling catas- 
trop he at the c luster core predicted by fluid mod- 
els llFabianI 1 19941: iBinnevI l2003l: JPeterson fc FabianI 120061 : 



iParrish. Quataert fc Sharmal2009l : lBogdanovic et al.H2009t ) 
This requires modelling various heating processes involving 
conversion of the energy of plasma motions (turbulent or 
otherwise) in to heat vi a some form of effective viscosity (e.g., 
Omma et al. 1[2004 ; Fa bian et al.ll2005bl:l iDennis fc ChandranI 



200^ 'Chandran fc Raseral '2007; Guo. Oh fc Ruszkowski 



Briiggen fc Scannapiecg ,2009. : iKunz et al.ll201ll ). the 
dynamical effect of thermal i nstabilities arising in the 



magneti zed ICM (lBalbuj|2000l; IParrish. Stone fc Lemasted 
2008; QuataertI l2008l: ISharma. Quataert fc Stond | 2008| : 



Sharma et all l2009l: IParrish. Quataert fc Sharma 



20ld: iBogdanovic et al .1 |2009| ; iRuszkowski fc Ohi boiC 
Schekochihin et al.ll201ol ). and the effective thermal conduc- 



200C 



tivity of this medium with acco unt taken of the tangled 
magnetic field |"Chandran fc Cowley.. 1998.: ,Malvshkin.2001 .: 
Naravan fc MedyedevI l200ll: IZakamska fc NaravanI booj ; 



Cho et al.ll2003l ; IVoigt fc FabianlbOol T 

• Can we construct theoretical and numerical models 
of the ICM dynamics that reproduce quantitatively the 
features we observe, e . g., the rise of radio bubbles (e.g., 
IRuszkowski et al.|[2007l : iDong fc Stonell2009l ). the formation 
and propagation of shocks, fronts and sound waves, the 
structure of ICM velocity, density, temperature fluctuations? 



In section [T] wc will discuss some of the relevant questions for 
astrophysical contexts other than galaxy clusters. In section [7.31 
wc will also give a brief survey of what in our view is the current 
state of play in answering the questions raised here in view of 
what we know about the plasma instabilities in the ICM and 
their likely saturation mechanisms. 



• Can we explain the origin of the 
netic fields ( probably via some form 



dynamo, see ISub ramanian. Shukurov 



cluster mag- 
of turbulent 



Haugcn '2006|: 

Schekochihin fc Cowley. 2006. : Enfilin fc Vogt.,2006i ; ,Xu ot all 
2009t ) and their observed spatial structure? 



Addressing these questions requires a theoretically 
sound mean-field theory for the ICM dynamics, i.e., a set 
of prescriptions for its effective transport properties (viscos- 
ity, thermal conductivity), which depend on the unresolved 
microphysics. Without such a theory, all we have is numeri- 
cal simulations based on fluid models (see references above), 
which, while they can often be tuned to produce results that 
are visually similar to what is observed, are not entirely 
satisfactory because they lack a solid plasma-physical basis 
and because reflning the numerical resolution often breaks 
the agreement with observations and requires retuning. A 
satisfactory transport theory is lacking because any plasma 
motions in the ICM that change the strength of the magnetic 
fleld trigger microscale plasma instabilities (see sections 12.21 
and l2.3p and we do not know what happens next. 

How some of these instabilities arise and evolve is dis- 
cussed in greater detail below. In order to make this discus- 
sion more quantitative, we need to flx a few physical param- 
eters that characterize the ICM. In reality, these parameters 
vary considerably both between different clusters and within 
any individual cluster (as a function of radius: from the 
cooler, denser core to the hotter, more diffuse outer regions). 
However, for the purposes of this discussion, it is sufficient to 
adopt a set of fiducial values. Let us consider the plasma in 
the core of the Hydra A cluster (also used as a representative 
example in our precedin g pap ers, ISchekochihin fc CowlevI 
I2OO6I: ISche kochihin ct al.' 2008), where the parameters are 
(|David et al. 2001; Enfilin & Vogt 200ii) 



• particle (ion and electron) number density 



rii = rie ~ 6 X 10 ^ cm 



• measured electron temperature is 
Te ~ 3 X 10^ K; 



(1) 



(2) 



the ion temperature is unknown, but assumed to be compa- 
rable, Ti ^ Te\ then the ion thermal speed is 

vtM = ( — 1 ~ 7 X 10^ cm s"' (3) 



A 



rrii is the ion mass, Ti is in erg); the ion Debye length is 



Vth: 



. 2 \ -1/2 



2 X 10'' cm; 



(4) 



• the ion-ion collision frequency (in seconds, assuming rii 
in cm~^ and Ti in K) is 



u^^ ~ 1.5n,r, ~ 5 X 10~" s-^ 



consequently the mean free path is 

mfp 



, Vt\^i 20 

AmfD ~ ~ 1.3 X 10 cm 



(5) 



(6) 



the rms magnetic field strength is (|Vogt fc Enfilin I 

(7) 



l2005h 

B ~ 7 X 10"'' G; 
consequently the plasma (ion) beta is 
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A = ^ ~ 130, 

the ion cyclotron frequency is 

a = — ~ 0.07 

rmc 



(8) 



(9) 



(e is the elementary charge, c the speed of light) and the ion 
Larmor radius is 

= ~ 10« cm; (10) 

note that the magnetized-plasma condition pi <C A^fp is 
satisfied extremely well; 

• the typical velocity of the plasma motions is 

[/ ~ 2.5 X 10^ cm s'^ (11) 

(cf. ISanders et al.l l2010bl . 1201 ll . who consider a sample of 
clusters), while the typical length scale of these motions is 



i ~ 2 X 10^^ cm; 

consequently the Mach number is 

M = — ~ 0.3 



(12) 



(13) 



(so the motions are subsonic, hence approximately incom- 
pressible on scales smaller than that of the mean density 
variation) and the Reynolds number based on coUisional par- 
allel viscosity is 



Re = 



LU 



Amfpt'thi 



60, 



(14) 



assuming Kolmogorov scalings for turbulence, the viscous 
cutoff scale is 



I ~ LRe"^/-* ~ 10^^ cm 

and the typical velocity at this scale is 

« ~ f/Re"^/* ~ 10^ cm s'V 



(15) 



(16) 



so the approximate rms rate of strain (assuming a viscous 
cutoff for the motions) is 



70 ~ y ~ R-e ~ 10 s . 



(17) 



2.2 Origin of pressure anisotropy 

If we consider length scales greater than p; and time scales 
longer than Q.i (which is easily true for any large-scale dy- 
namical processes in the ICM), the momentum equation for 
the plasma flow, characterized by the mean velocity u, is 
(e.g.,[K^srud, 198j) 



du 

'd^ 



V 
V 



Stt 



bb 



-P\\ 



(18) 



where du/At = i9/9t + u ■ V is the convective derivative, 
b is the unit vector in the direction of the local magnetic 
field, B is the field's strength, and p±_ and are the per- 
pendicular and parallel plasma pressure, which are the only 
components of the plasma pressure tensor that survive at 



these long spatial and temporal scales: 



rus j A'vvv fs = p^[\-bb'j +p^\bb, (19) 

d^v^/., (20) 



j3 2 J. 



(21) 



where fs is the distribution function for species s {s = i,e), 
V its velocity variable (particle's peculiar velocity), and v± 
and W|| the projections of v perpendicular and parallel to 
the magnetic field. The magnetic field is determined by the 
combination of Faraday's and Ohm's laws, which at these 
long scales takes the form of the ideal induction equation 



dS 

dt 



B Vu BV u. 



(22) 



Without as yet going into the technicalities of kinetic 
theory, it is not hard to show that pressure anisotropics arise 
naturally in a weakly collisional plasma. Indeed, the first 
adiabatic invariant p = v^/2B of a gyrating particle is con- 
served on time scales intermediate between the collision time 
and the cyclotron period (a nonempty interval when plasma 
is magnetized, i/u <^ Hi). Since p± is proportional to the 
sum of the values of jj, for all particles, p±/B should be a 
conserved quantity, i.e., if the magnetic field changes (as a 
result of plasma motions into which the fiux is frozen, see 
equation (|22|l 'l then p± should change accordingly. For the 
purposes of this qualitative discussion, we may momentarily 
ignore the fact that changing B also causes to change 
(in a different way from p± ; see Appendix IA2.15[l and so 
conclude that changing B will cause pressure anisotropics 
to develop. 

In the absence of collisions, the pressure anisotropics 
would track the field strength. If collisions do occur, 
even weakly, their effect will be to relax the system to- 
wards an isotropic pressure (and a Maxwellian distribution). 
Thus, there is a competition between changing B inducing 
anisotropy and collisions causing isotropization. This can be 
modelled by the following heuristic equation: 

1 dpx 1 dB p± — P|| 
p± dt B dt " p± 



bb : Vu 



-P\\ 



Pi- 



rn) 



where we have used equation p2|l to express the change in 
the field strength in terms of the plasma fiow velocity and 
assumed, for the purposes of this qualitative discussion, that 
plasma density is constant (i.e., the motions are incompress- 
ible). Considering what happens on time scales longer than 
the collision time, we conclude, after examining the right- 
hand side of equation p3|l , that we should expect the typical 
(ion) pressure anisotropy in a moving plasma to be 



■P\\ 



1 1 dB 

Vii B dt 



To 



(24) 



where 70 is the typical rate of strain of the plasma motion 
Thus, the pressure anisotropy is regulated by the ratio of 



A few tangential comments are appropriate here; 
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the typical rate of change of the magnetic-field strength to 
the collision frequency. 

Substituting the numbers from section [2m we find that 
I A| ~ 0.02 in the core of Hydra A. Is this a large number? It 
turns out that it is a huge number because such anisotropies 
will make the plasma motion violently unstable. 



2.3 Firehose instability 

While the full description of the plasma instabilities trig- 
gered by pressure anisotropies requires kinetic treatment, it 
is extremely straightforward to deduce the presence of the 
firehose instability directly from equation (|18|l . 

Consider some "fluid" solution (uo, Bo,Po±,Po||) of 
equations (|18[l and (|22[) that varies on long time and spatial 
scales — that can be thought of as the turbulence and/or 
some regular magnetofluid motion caused by global dynam- 
ics. Let us now examine the linear stability of this solution 
with respect to high-frequency {u S> lVno|), short-scale 
(fe S> |Vno|/tio) perturbations (Su, SB , Sp±, Sp^^). Mathe- 
matically, this is simply equivalent to perturbing a straight- 
magnetic-field equilibrium of equations (|18p and (|22[) : 



miUiuSu = — fcx I 5p 



4-K 



+ fc||56 ( Pox - Poll + 



op II + (po± - Poll) 



k^^Su — bo (k ■ Su) 



(25) 
(26) 



where 5b = SB±/Bo, we have used k ■ 5b — — fc||i5i3||/_Bo 
(from V • i? = 0), and _L and || are with respect to the 
unperturbed magnetic field direction bo. Pressure perturba- 
tions can only be calculated from the line arized kinetic equa- 
tion (see, e.g.. ISchekochihin et al.ll2005l ). but even without 



(i) The electron pressure anisotropy is smaller by a factor of 
~ 43 because the electron collision frequency is ~ {mi/me)^^'^L'ii. 

(ii) If we use equation I I23II to write explicitly p± — py = 
{pj_/uii)bb : Vii and substitute this into equ ation II18II. we re - 
covcr (to lowest order in ua/fli) the well known iBraginskiil lll965l ) 
momentum equation with anisotropic viscosity, where p±/uii ~ 
miriiV^-^^/vii is the Braginskii parallel viscosity coefficient. 

(iii) If a Kolmogorov-style turbulence is assumed to exist in 
the ICM, the typical rate of strain 70 will be dominated by the 
motions at the viscous cutoff scale. However, as we saw in sec- 
tion l2.ll the Reynolds-number estimates for ICM do not give very 
large values and one might wonder whethe r calling these motions 
turbulence is justified jPabian et al .l2003bh . However, for our pur- 
poses, it is not important whether the rate of strain is provided 
by the viscous cutoff of a turbulent cascade or by a single-scale 
motion because either ca n change the magnetic field a nd thus 
cause pressure anisotropy jSchekochihin fc Cowlevll200(J) . 

(iv) For a purely compressive motion, A ~ —V ■ u/Zva (i.e., 
the anisotropy is still related to the change in the magnetic-field 
strength; see equation I I22II ) , but one has to work a little harder to 
show this. In the compressible case, one also discovers that heat 
fluxes contribute to the anisotropy alongside velocity gradients 
(this is done in Appendix I A2.131 sec equation IIA63I1 ). 



knowing them, we find that for the Alfvenically polarized 
modes, 5ti oc bo x k, the dispersion relation is 



, , , POX - Poll , 2 

mini 



1/2 



/3 



1/2 



(27) 



(pox/miTii)^''^, A = (po 



where va = Bo /\/4:iTmini, Ca 
Po||)/Po± and /3 = 87rpo±/Bo- 

Equation (|27|l is simply the dispersion relation for 
Alfven waves with a phase speed modified by the pres- 
sure anisotropy. If the pressure anisotropy is negative, 
A < 0, the associated stress opposes the Maxwell stress 
(the magnetic tension force), the magnetic-field lines be- 
come more easily deformable, the Alfven wave slows down 
and, for A < —2/13, turns into a nonpro pagating unstable 



mode — this is the firehose instabi li tv (iRosenbluthI 
Chandrasekhar. Kaufman fc Watson 19581: Parked 



Vedenov fc Sagdeev 
196lh . Its growth rate can. 



1956 



1958 



19581 : IVedenov. Velikhov fc SagdeevI 
in general, be almost as large as 
the ion cyclotron frequency as fcypi approaches finite values 



(see section \2.b\ . For the ICM parameters given in section 
2.11 the instability is, therefore, many orders of magnitude 
faster than either the large-scale dynamics (typical turnover 
rate ~ |Vuo| ~ 70) or collisions (typical rate uu). 

Thus, any large-scale motion that leads to a local de- 
crease in the strength of the magnetic fielcQ gives rise to 
a negative pressure anisotropy, which, in turn triggers the 
firehose instability, producing Alfvenically polarized fluctu- 
ations at small parallel scales — unless the plasma beta is 
sufficiently low (magnetic field is sufficiently strong) for the 
magnetic tension to stabilize these fiuctuations. Using the 
typical size of A estimated at the end of section [ST^] for the 
Hydra A ICM parameters, we find that the typical beta be- 
low which the firehose is stable is /3 ~ 100, which is quite 
close to the measured value (see section [TT} — perhaps not 
a coincidence? 

Positive pressure anisotropi es also lead t o instabilities 
most importantly, mirror; see iFurth' 'l9 6j: iBarnes} 
■ 1; [SoiTtl 



Taiiril Il967l: iHasegawal I l969l : ISouthwood fc KivelsonI 



1966 



1993 



HeUingeill2007[ and references therein), but they involve res- 



onant particles and are mathema tically harder to handle . 
We will not discuss them here (see ISchekochihin et al.ll2005l . 
I2OO8I : iRincon. Schekochihin fc Cowlevll2010l) . 



2.4 Nonlinear evolution of the firehose instability 

A nonlinear theory of the firehose instability can be con- 
structed via a quasilinear approach, in which the unsta- 
ble small-scale (perpendicular) fluctuations of the mag- 
netic field on the average change the local magnetic-field 
strength and effectively cancel the pressure anisotropy 
jSchekochihin et al]|2008l ') . In equation p4)) , let us treat the 
changing magnetic field as the sum of the large-scale field 



* While turbulence on the average is expected to lead to the 
growth of the magnetic fie ld (the dynamo effect; see, e.g., 
ISchekochihin fc Cowlevl2006l and references therein) , locally there 
will always be regions where the field strength (temporarily) de- 
creases. Decrease of the field and, consequently, negative pres- 
sure anisotropy can also result from expanding motion, which 
decreases the density of the plasma — as, e.g., in the solar wind. 
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and the small-scale firehose fluctuations: B — Bo + SBj 
Then the field strength averaged over small scales is 



B ~ Bo 1 + 



1 \5Bj 



^0 



(28) 



where the overbar denotes the average (under which small- 
scale fluctuations vanish). The contribution from SB± is 
small, but for large enough fcy, it is growing at a greater 
rate than the rate of change of the large-scale field, so its 
time derivative can be comparable to the time derivative 
of Bo- As Bo is assumed to be decreasing, the growth of 
the fiuctuations can then cancel this decrease and drive the 
total average pressure anisotropy to the marginal level, A = 
—2/13. From equation (|24p . we get 



v^^ \Bodt 2 At 



d \5Ba 



2 



(29) 



The rate of change of Bo is the typical rate of strain of the 
(large-scale) motion, (l/Bo)dBo/df ~ — l7o|- The firehose 
growth rate 7 = —iuj is given by equation (|27p . As long as 
the firehose fiuctuations are smaller than the critical level 



\5B^ 



-"0 



1 70 1 

7 



(30) 



they cannot enforce the marginality condition expressed by 
equation (|29p and will continue growing until they reach the 
required strength (which is still small compared to the large- 
scale field because I70I/7 ^ 1 for sufficiently large fcy). After 
that, their evolution becomes nonlinear and is determined 
by equation (|29|) . whence we find that their energy has to 
grow secularly: 



\5Ba 



^0 



70 



t. 



(31) 



As long as the large-scale field keeps decreasing, the small- 
scale fiuctuation energy cannot saturate because if it did, 
its time derivative would vanish, the anisotropy would drop 
below marginal and the instability would come back. 

The secular growth given by equation (|3ip leads to 
5B±/Bo ~ 1 after roughly one turnover time (~ |7o|~^) 
of the large-scale background motion that produces the 
anisotropy in the first place — thus, the magnetic field can 
develop order-unity fiuctuations before this background mo- 
tion decorrelates. What all this means for the large-scale 
dynamics on longer timescales, we do not know. 

In what follows, we will be guided by the simple ideas 
outlined above in constructing a more rigorous kinetic the- 
ory of the nonlinear firehose instability. 



in which equation (|18|) is derivedlf] so finite-Larmor-radius 
(FLR) corrections must be brought in. 

Direct calculation of the linear firehose growth rate 
from the hot-plasma dispersion relation shows that the 
peak of the growth rate is at fcypi ~ |A + 2/(3\^^^ 
for the parallel {k± = 0) firehose (|Kennel fc SagdeevI 
1 19671 : iDavidson fc Volkl Il968l . this result will emerge in 
section 14. 2p and, in genera l, at kpj ~ 1 for the 
oblique firehose with k± ^ (|Yoon. Wu fc de Assislll99j : 
iHellineer fc Matsumotol I2000I ). This means that the maxi- 
mum growth rate of the instability is 7max ~ | A -f 2//?|Sli ~ 
10^'^ s^^ for k± — (see section 14. 2p and 7max ~ |A + 
2//3|^/^ni ~ 10"^ s"^ for fcx / 0, where we have used the 



ICM parameters of section 12.11 and the estimate of A from 
section 12.21 

There are two conclusions to be drawn from this. First, 
the linear instability is enormously fast compared with the 
large-scale dynamics that cause it, so its nonlinear behaviour 
must be fundamentally important at all times. Second, in or- 
der to understand the spatial structure of the firehose fluctu- 
ations, we need a theory that takes the FLR effects explicitly 
into account because it is the FLR that sets the scale and 
the growth rate of the fastest-growing mode. We now pro- 
ceed to construct such a theory for the simplest case — the 
parallel {k± = 0) firehose instability. 



3 KINETIC THEORY 
3.1 Basic equations 

The distribution function fs{t,r,v) satisfies the Vlasov- 
Landau kinetic equation 



Ot rUs \ c 



dv 



C[fs 



(32) 



where s = z, e is the particle species, r its position, v ve- 
locity, Qs and nis are the charge and mass of the particle of 
species s (qe = — e, qi = Ze, Z = 1 for hydrogen plasma), 
E and B are the electric and magnetic fields, and the term 
on the right-hand side is the collision operator. The electric 
and magnetic fields are determined from Maxwell's equa- 
tions: quasineutrality 

qsTis = "^18 I d^i; /s = (33) 

s s ■' 

{us is particle number density). Ampere's law 

j = ^ qsUsUs = '^<ls d^v V fs = ^ V X B (34) 



2.5 Effect of finite Larmor radius 

We have so far carefully avoided discussing the magnitude of 
the wavenumber fcy of the firehose fiuctuations, simply refer- 
ring to them as "small-scale," with the implication that their 
scale would be smaller than that of the background fiuid dy- 
namics that cause the instability. Examining the dispersion 
relation (|27p . we see that the growth rate of the instabil- 
ity is proportional to fcy , so the smaller the scale the faster 
the instability. This ultraviolet catastrophe cannot be re- 
solved within the long-wavelength approximation, kpi <§; 1, 



{j is current density, Ug is the mean velocity of the species 
s), Faraday's law 



dB 

dt 



-cV X E, 



(35) 



and V ■ B = 0. Note that equations (|33|l and (|34p are valid 
as long the particle motion is nonrelativistic and the scales 
we are interested in are larger than the Debye length. 



^ Which means that the equation is ill posed and eannot be solved 
without some kinetic prescription for the handling of small scales. 
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It is convenient for what follows to calculate the dis- 
tribution function in terms of peculiar velocities v' = v — 
Us{t,r). Transforming the variables {t,r,v) — >■ {t,r,v'), we 
find that equation (|32|l takes the form 



dt 



■ V/. + v' ■ Vfs 



+ 



dus 
'df 



^ Us X B v' X B 
E+— + 



^=C[f4. (36) 



We will henceforth drop the primes, v will be the pecu- 
liar velocity in all that follows. In this new formulation, the 
strategy for solving equations p3H36p is as follows. 



3.2 



Electron kinetics: Ohm's law and induction 
equation 



The electron kinetic equation can be expanded in the square 
root of the electron-ion mass ratio {me/mi)^^^ ~ 0.02, a 
natural small parameter for plasma. This expansion is car- 
ried out in Appendix lAll where we also explain what as- 
sumptions have to be made in order for it to be valid. The 
outcome of the mass-ratio expansion is that electrons are 
Maxwellian If] isothermal (Te = const), and the electric field 
can be determined in terms oiue, B and via a generalized 
Ohm's law: 

UeX B _ Vpe _ TeVn. 

erie 



E + 



c erie erie 

This can now be recast in terms of moments of the ion dis 
tribution: from equation (|33p . 

rie = Zm 



and from equation (f3^ 



3 

erie 



Ui 



47renf 

so equation (|37|l becomes 



V X B, 



E + 



u^x B _ TeVui ^{V X B)x B 

c erii AixZeni 



(38) 



(39) 



(40) 



and Faraday's law H35I) takes the form of the standard in- 
duction equation with a Hall term: 



dB 

dt 



= V X 



Ui — 



AnZerii 



V X B X B 



(41) 



3.3 Ion kinetics: continuity and momentum 
equations 

To close this set of equations, we must determine rii and 
Ui. Integrating equation (|36|) . we find that Ui = J di^v fi 
satisfies the continuity equation 



^ -f V ■ (n,nO = 0. 



(42) 



^ This means they do not contribute to the pressure anisotropy, 
which, to lowest order in the mass ratio, they indeed should not 
do, as pointed out already in footnote [3] Note that the valid- 
ity of these statements depends on the ordering of the collision 
frequencies given by equation ||A3|| . 



The equation for Ui (the ion momentum equation) follows 
from equation ((36]) for s = i, by taking the v moment and 
enforcing J d^v v fi — (by definition of the peculiar veloc- 
ity v), which gives 



dui 
'dt 



+ Ui ■ Vui 



V ■ ZTeVn, ^ {V xB)xB ^^g^ 



mini 



miTii 



Aumirii 



where we have used equation (|40p . the second term on the 
right-hand side is the electron pressure gradient, and we have 
introduced the ion pressure tensor 



Pi = mi I d V vv fi. 



(44) 



It is in order to calculate Pi in terms of Ui and B that we 
must solve the ion kinetic equation. We do this by means of 
an asymptotic expansion in a physical small parameter. 



3.4 Asymptotic ordering 

The small parameter we will use is expressed in terms of 
the M ach and Reynolds numbers (jSchekochihin et al.ll2005l . 

-5^71-0.1, (45) 

where we used the ICM parameters of section 12.11 This is 
the natural small parameter for the plasma motions because, 
using equations (|11H17[) . it is easy to see that 



(37) 



X 



mfp 



Vth; 



(46) 



where / is the viscous scale and u the typical flow velocity at 
this scale. The typical rate of strain 70 ~ u// is the relevant 
parameter for determining the size of the pressure anisotropy 
because, even though the viscous cutoff we are using is based 
on the parallel coUisional viscosity and so motions can exist 
below this scale, these mot ions do not change the streng th 
of the mag netic field fsee ISchekochihin fc CowlevI |2006| )PI 
Thus, the pressure anisotropy is (from equation (|24p ') 



70 
Vii 



u \ 



mfp 



(47) 



We solve the ion kinetic equation by asymptotic expan- 
sion in e. All ion quantities are expanded in e, so 



fi ~ foi + fli + f2i + f-ii + ' ' 

rii — noi + nii + n2i +713^ + - 

Ui = um + uii H , 

B = Ba + Bx+ . 



(48) 
(49) 
(50) 
(51) 



^ As already pointed out in footnote [S] our considerations do not 
depend on Re being large. If a single-scale flow is considered, our 
expansion is simply an expansion in Mach number. 
* This statement applies to macroscopic motions: for example, 
Alfvenic turbulence below the parallel viscous scale that can oc- 
cupy a wi de range of scales all the way down to the ion Larmor 
scale (e.g.. [Schekochihin et al]|2009l ). The fast, microscale plasma 
fluctuations triggered by plasma instabilities, including the fire- 
hose fluctuations that will be considered in this paper, will, on 
the average, change the field strength (see section [2.411 . Accord- 
ingly, their ordering [equation I I55II 1 will be arranged in precisely 
such a way that they are able to have an efl'ect comparable to the 
macroscale motions that produce 70. 
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The lowest-order quantities noi, uoi, Bo are associated with 
the motions that produce the pressure anisotropy and have 
the length scale I and time scale 70, so we order 



70 



(52) 



Since the instability parameter is A + 2//3i, we must order 
Bo so that 



2 

Ji 



Bo 



(53) 



The perturbations nn, uu, Bi around this slow large-scale 
dynamics are assumed to be excited by the prallel (fcj_ = 0) 
firehose instability and have much shorter spatial and time 
scales. Their typical wavenumber is the one at which the 
instability's growth rate peaks and their time scale is set by 
this maximum growth rate (see section [275)) : 



1/2 



(54) 



In order to be able to proceed, we must order the time 
scales of the lowest-order ("equilibrium") fields and of the 
fluctuations with respect to each other. Physically, they de- 
pend on different things and are not intrinsically related. 
However, our a priori consideration of the nonlinear evo- 
lution of the instability (section 12. 4|) suggests that for the 
nonlinearity to become important, we must have (see equa- 
tion ((soil) 



Bi 

Bo 



1/2 



Since Bi ~ eBo, this tells us that we must order 

2 4o ^ 2o 2\ 

70 ~ e 7 ~ e i2i => Vii e ih, pi ^ e Xmfp. 



(55) 



(56) 



These relations are, of course, not strictly right in the quan- 
titative sense — the Larmor radius is grossly overestimated 
here if we take the value of e for the ICM given by equation 
H45p and then compare what equation (|56|) gives us as the 
value of pi with the ICM estimate in section 12.11 (equation 
(|10|) '). However, ordering pi this way allows us to capture all 
the important physics in our formal expansion. We will also 
argue in section [4.3.21 that this ordering of the finite Larmor 
radius physics gets quantitatively better as the nonlinear 
regime proceeds (see footnote I13p . By the same token, the 
growth rate of the instability in the ICM is typically much 
larger than the collision rate, while we have ordered them 
similar — but again, this ordering formally allows all the 
important physical effects to enter on a par with each other 
and also gets better in the nonlinear regime, where the fire- 
hose fluctuations grow slower. 

Let us summarize our ordering of the relevant time and 
spatial scales compared to fcyVthi and fcy, respectively: using 
equations H56I) and (|54|) . we have 



70 ~ e fcyUtiM, 



Vii ~ efcyVthi 



mfp 



Pi 



(57) 
(58) 



3.5 Firehose fluctuations 



The ordering we adopted, inasmuch as it concerns the prop- 
erties of the firehose fluctuations, applies to the parallel flre- 
hose only, so we now explicitly restrict our consideration to 



the case of Vx = for all flrst-order perturbations. Since 
V ■ B = 0, this immediately implies 



BI = 0, 



(59) 



so Bi = . Here and in what follows, || and _L refer to 
directions with respect to the unperturbed field Bq. 

The induction equation (|4ip. taken to the lowest order 
in e, gives 



(60) 



(all terms here are order e^fc||i;thi; see section [3741 note that 
the Hall term in equation (1411) is subdominant by two orders 
of e). Here A /At — d/dt + uo ■ V is the convective derivative, 
but, since Vuo ~ e'^fc||?;thi, the shearing of the perturbed 
field due to the variation of uq is negligible and we can 
replace d/dt by by transforming into the frame moving 
with velocity uq. 

In the continuity equation (|42p taken to the lowest order 
in e, setting Vx = gives 



d nu II 
at noi 



(61) 



(all terms are order e^k^^vthi)- Anticipating the form of the 
unstable perturbation, we will set 

nu = 0, ul, = (62) 

without loss of generality. In Appendix lA2.6l we will explic- 
itly prove that nu = 0. In Appendix lA2.9l we will learn that 
n2i = as well. 

Consider now the ion momentum equation (|43p . In the 
lowest order of the e expansion (terms of order efcuWthi)! it 
gives, upon using equation 



V Pu = 0. (63) 

We will learn in Appendix IA2.8I that this can be strength- 
ened to set 

Pi. = 0. (64) 

In the next order (e^k^^Vthi), we get (using n2i = 0) 

V-Po, + V-P2» + ^reVno, = 0. (65) 

Averaging this over small scales eliminates the perturbed 
quantities, so we learijfl 



V ■ Poi + ZnVno, = 

and, therefore, from equation (|65p . also 

V ■ P2, = 



(66) 



(67) 



(confirmed in Appendix IA2.9p . Finally, in the third order 
(e^fc||i;tijj), the perpendicular part of equation (|43p deter- 
mines the perturbed velocity field: 



dui, (y_Psi)^,2^ Bi 

dt niinoi Bo 



(68) 



^ This is simply tlic pressure balance for the large-scale dynamics, 
an expected outcome for a system with low Mach number. In 
Appendix IA2.5I we will show that the zeroth-order distribution 
is Maxwellian, so the pressure associated with it is a scalar, poi = 
noiToi, and equation I I66I I becomes (To; +ZTe) Vnoi -|-noi VTo; = 
0. Further discussion of the role played by the ion temperature 
gradient can be found in section [5] 
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where va ~ Bq /Anminoi. There is no ZTe^^n^i term in 
equation (|68|l because we assume that the only small-scale 
spatial variations of all quantities are in the parallel direc- 
tion. The ion pressure term (V • PaOx calculated 
by solving the ion kinetic equation (see section [377)) . 

To summarize, we are looking for perturbations such 
that Vx — 0, nii — 0, = 0, u^^ = 0, while Bj^ and 
uj^i satisfy equations (|60|l and (|68|) . Physically, this reflects 
the fact that the parallel {k± — 0) firehose perturbations 
are Alfvenic in nature (have no compressive part). That it 
is legitimate to consider such perturbations separately from 
other types of perturbations is not a priori obvious, but will 
be verified by our ability to obtain a self-consistent solution 
of the ion kinetic equation, which will satisfy equations (|64p . 
((Sni), and ((571) (see Appendix IX2| . 

3.6 Large-scale dynamics 



In section 13.51 equations for the first-order fields, Un and 
Bi emerged after expanding the induction equation (jdip 
and the continuity equation (|42p to lowest order in e and 
the momentum equation (|43p up to the third order. If, using 
the ordering of sect ion [HTil we go to the next order and aver- 
age over small scales to eliminate small-scale perturbations, 
we recover the equations for the large-scale (unperturbed) 
fields: the induction equation 



dt 



(all terms are order e'^k^^VthiBg), the continuity equation 



dt 



-noiV ■ uoi 



(70) 



(all terms are order e k^ivthinoi), and the momentum equa- 
tion 

duoi _ _Bi Bo VBo 

m-iTioi— t;- = -V ■ P2— V- \ (71) 

dt Stt 47r 

(all terms are order e'^rriinoik^^v'^Y^i). The divergence of the 
second-order ion pressure tensor here is with respect to the 
large-scale spatial variation (according to equation (|67p . it 
has no small-scale dependence). Again, P2i is calculated 
from ion kinetics. 

Equations (|69H7ip are precisely the kind of mean-field 
equations that are needed to calculate the large-scale dy- 
namics of astrophysical plasmas. They look just like the 
usual fluid MHD equations, the only nontrivial element be- 
ing the pressure term in the momentum equation (|71|l . The 
goal of kinetic theory is to calculate this pressure, which 
depends on the microphysical fluctuations at small scales. 
In this paper, we only do this for the parallel {k± = 0) 
fir ehose fluctuations. For the mirro r fluctuations, it is done 
in iRincon. Schekochihin fc CowlevI (|2010D (using a some- 
what different, near-marginal-stability asymptotic expan- 
sion), while the oblique flrehose fluctuations are a matter 
for future work. The implications of our results for the ion 
momentum transport will be discussed in section 14.41 

3.7 Solution of the ion kinetic equation 

We now proceed to use the ordering established in section lT^ 
to construct an asymptotic expansion of the ion kinetic equa- 
tion. This procedure, while analytically straightforward, is 



fairly cumbersome and so its detailed exposition is exiled to 
Appendix I A2I The results are as follows. 

In the expansion of the ion distribution function (equa- 
tion (|48|) '). foi is found to be a Maxwellian f Appendix lA2.5|) . 
with density noi and temperature Toi that have to satisfy the 
equilibrium pressure balance constraint (see equation (|66|l 
and Appendix IA2.8|) . 

The flrst-order perturbed distribution function, fu, is 
proportional to bo- VToi and is responsible for the large-scale 
coUisional ion heat fluxes (Appendix IA2.8|I . 

The second-order perturbed distribution function f2i 
contains the pressure anisotropy. The corresponding second- 
order pressure tensor is diagonal: 



P2» =P2i\ + (P2» 



(72) 



where p2i is the perturbed isotropic pressure and P2i — P^i 
the lowest-order pressure anisotropy. The isotropic part of 
the pressure is determined from the large-scale equations for 
density noi, temperature Toi and velocity uoi of the fluid — 
this is explained in detail in Appendix lA2.12l but here let us 
just assume for simplicity that the zeroth-order density and 
temperature are constant (Vnoi = 0, VToi = 0), in which 
case the continuity equation (|70|l reduces to V • uoi = 
and p2i then follows from enforcing this incompressibility 
constraint on the momentum equation (|71|l . The pressure 
anisotropy is calculated in Appendix IA2.131 



A(t) 



■Pii 



POi 



Ao + 



■ f At'e 
Jo 



dt' 



^0 



(73) 



where poi = noiToi is the equilibrium pressure, the overbar 
denotes the averaging over small scales of the nonlinear feed- 
back on the anisotropy from the firehose fluctuations, and 
Ao is the pressure anisotropy arising from the large-scale 
motions. In general, it contains contributions from changes 
in the magnetic field strength (because of the approximate 
conservation of the first adiabatic invariant, as discussed 
qualitatively in section I2.2|l . compression and heat fluxes 
(see equation ()A63|) ). When noi and To; are constant, only 
the anisotropy induced by the changes in fleld strength sur- 
vives which is the case we will consider here: 



Ao = 



bubo : 'Vupi 

Vii 



1 1 dBo 

Uii Bo dt 



To 



(74) 



This is exactly what was anticipated qualitatively — see 
equation (|24p . Since we are interested in the flrehose insta- 
bility, we assume Ao < oF^ 

Finally, the third-order perturbed distribution function 

^'^ The effect of heat fluxes on the firehose turbulence is consid- 
ered in section [5] 

We have assumed the initial anisotropy A(0) = Ao. It is 
equally possible to start from any other value, including stable 
situations. In that case, the large-scale drivers of the pressure 
anisotropy will gradually build it up to the maximum (nega- 
tive) level, Ao, whereupon further evolution will proceed in the 
same way as discussed below. Mathematically, this amounts to 
mutiplying Ao in equation II73I I by (1 — e"''""') (solution of 
equation IIA6OI1). An exa mple of such a set up can be found in 
ISchekochihin et al.l 1 I2OO8I) . Note that if the initial fluctuation level 
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/si is responsible for the third-order pressure tensor that 
appears in the perturbed ion momentum equation (|68|) . 
The relevant part of that tensor is calculated in Appendix 
IA2.14I Assuming constant density and temperature (other- 
wise, there is again a contribution from the heat fluxes; see 
Appendix [5|, it may be written as follows 



(V- Pa 



= -V|i 



POi 



-DO Hi 



(75) 



where A{t) is given by equation (|73|l . 

Let us now use these results to study the firehose tur- 
bulence (sections I4.1H4.3P and its effect on the large-scale 
dynamics (section [J^J . 



4 FIREHOSE TURBULENCE 



4.1 Firehose turbulence equation 



Using the results derived in Appendix I A2I and summarized 
in section ISTfl we find that the ion momentum equation (|68|) . 
which describes the evolution of the perturbed ion velocity, 
is, in the reference frame moving with uoi, 



dt 



Vthi 



A{t) + 



/3. 



-Bo 



(76) 



where we have used equation (|75|) for the pressure term in 
equation ((68}. Three forces appear on the right-hand side of 
this equation. First, there is the stress due to the anisotropy 
A of the ion distribution, given by equation (|73|) . The lat- 
ter equation is the quantitative form of the expression for 
A(t) that we guessed in equation (|29|) : the first term in 
equation (|73|) is due to the slow decrease of the large-scale 
magnetic field, the second to the average effect of the grow- 
ing small-scale fiuctuations, which strive to cancel that de- 
crease. The second term in equation (|76|) . proportional to 
l//3i = v\/v^^^i, is the magnetic tension force, which resists 
the perturbation of the magnetic-field lines and, therefore, 
acts against the pressure-anisotropy driven instability. The 
instability is marginal when A -I- 2//3i — >■ —0. Finally, the 
third term is the FLR effect, which, as was promised in sec- 
tion [2]5] and as will shortly be demonstrated, sets the scale 
of the most unstable perturbations. 

Let us now combine equation (|76|l with the induction 
equation (|60p for the perturbed magnetic field, also taken 
in the reference frame moving with uq. After differentiating 
equation H60[) once with respect to time, we get 



1 dBj 

fii dt 



X be 



(77) 



In the second term on the right-hand side, we have used 
equation (|60|l to express W\\Un in terms of the time deriva- 
tive of Bi . Equation (|77|l with A(t) defined by equation 
H73p is a closed equation for the perturbed magnetic field 
with nonlinear feedback (last term in equation ((73])). This is 



is not infinitesimal, the nonlinear quenching of the anisotropy 
(discussed in the subsequent sections) can start before the maxi- 
mum anisotropy is built up. 




Figure 1. Proqucncies (thin lines) and growth rates (bold 
lines) of the unstable firehose modes (red/solid: the "+" mode; 
blue/dashed: the "— " mode) given by equation jTOj. The insta- 
bihty parameter here is A -|- 2//?^ = —0.01. Dotted vertical lines 
indicate the wavenumber of fastest growth kp = 0.2 (equation 
(|8TJ) and the dotted horizontal lines the corresponding maximum 
growth rate 7max = ImWp = 0.01 (equation Il82|l '). 



the equation for the one-dimensional (fcx = 0) firehose tur- 
bulence. It represents the simplest nonlinear model for this 
kind of turbulence available to date0 



4.2 Linear theory 

In the linear regime, we may neglect the second term in 
equation ((73}, so A = Aq. The linear dispersion relation for 
equation (I77p is 



2 ^il"th« / ^ 2 

uj — - — Ah — 

2 V A 



(78) 



This has four roots out of which two are unstable when 
A + 2/ Pi < 0: 



UJ . . \k\ 

n> 4 72 



A + 



1/2 



fe2 



where k — k^^pi and 

1/2 



ko = 2V2 



(79) 



(80) 



(this linear dispersion relation was first ob tained by 
iKennel fc Sagdeevll 19671 : (Davidson fc Volklll968l ). Unlike in 



The essential differen ce with the equation we derived in 
ISchekochihin et al.l 1I2OO8I ) is the FLR term, which removes the ul- 
traviolet catastrophe of the long- wavelength firehose and thus al- 
lows equation l|77ll to handle non-monochromatic (multiscale) so- 
lutions. In scction r4.3l we will see that this produces a mu c h mor e 
complex behaviour than was seen in ISchekochihin et al.l l|2008h . 
justifying the term "firehose turbulence." 
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Figure 2. Left panel: evolution of the magnetic energy \Bj^\'^ / Bq = ^Zfe l^feP with time in a numerical solution of equations l|86|l and 
II87I I with parameters 1981 1 ; the time here is normalized using the collision frequency va , not the cyclotron frequency ; the two horizontal 
lines show the "coUisional" (lower line) and "coUisionless" (upper line) estimates for the energy at which the nonlinear feedback turns on: 
equations II89I I and 1901 1 , respectively; the red dotted line shows the nonlinear asymptotic given by equation 1911 1 . Right panel: evolution of 
the instability parameter (pressure anisotropy) A + 2//3i in the same numerical solution. Inset: log- log plot of the evolution of |A + 2//3i|; 
the red line shows the slope corresponding to 1/t (see equation 195l l). 



the long- wavelength limit (fcypi — 0), there is now a real 
frequency (so the firehose perturbation propagates while 
its amplitude grows exponentially and the vector ro- 
tates; see section l4.3.ip and the growth rate has its peak at 



ko/V2, 



1/2 



= i±l + i) 



-I 



(81) 
(82) 



where the complex peak frequency ojp is in units of fii. At 
k\\Pi > ko, there is no growth and the firehose perturba- 
tions turn into purely propagating Alfven waves (modified 
by pressure anisotropy and dispersive FLR corrections). 

The dependence of the frequencies and growth rates of 
the two unstable modes on wavenumber given by equation 
(|79p is plotted in figure [T] for a representative value of the 
instability parameter A + 2/l3i — —0.01 (this is the value 
used in the numerical solution of section [4. 3. 3[l . 

It should be pointed out here that in this theory, there 
is no dissipation of the magnetic fluctuations excited by the 
firehose. The most unstable wavenumber is set by dispersive 
effects; the stable modes are undamped. 



4.3 Nonlinear evolution and spectrum 

4-3.1 Firehose turbulence equation in scalar form 

Since the nonlinearity involves the spatially averaged per- 
turbed magnetic energy, the firehose turbulence is compactly 
described in Fourier space not just in the linear but also in 
the nonlinear regime: this amounts to replacing Vy — > — fcy 

in equation dTT)) and \Bi\^ = Y^k^^ l-Si"(fc||)P in equation 
(|73p . A simple ansatz can now be used to convert equation 
(|77p into scalar form. Let 



A f''\ 



(83) 



= Ak(t)sm{—t + 4,k 



(84) 



where the axes {x,y) in the plane perpendicular to 6o 
are chosen arbitrarily and we have non-dimensionalized 
wavenumbers and time: 



k\\pi 



t. 



(85) 



This ansatz amounts to factoring out the rotation of 
the vector Bj^(k) (the first term in equation (|79p ). The 
wavenumber-dependent but time-independent phase (pk is 
determined by the initial condition. We assume (f>k = (l>-k, 
so Al = A-k must be satisfied to respect the fact that Bj^ 
is a real field. The fluctuation amplitude Ak{t) satisfies 



Ak, 



A{t) = Ao + - 



(86) 
(87) 



where = Vu/Qi = Pi/Amfp and we remind the reader 
that Ao < 0. It is manifest in the form of equation (|86|l 
how the dispersion relation (|79|l (without the first term) is 
recovered. Note that there is no coupling between different 
wavenumbers modes in the sense that if a mode is not ini- 
tially excited, it is never excited. The only effect that modes 
have on each other is via the sum over k in equation (|87|) . 
to which they all contribute. 



4-3.2 Qualitative picture 

Already on the basis of linear theory and the qualitative 
considerations of section we can construct a fairly clear 
picture of the evolution of the firehose turbulence. Assuming 
a broad-band infinitesimal initial perturbation in k space, 
at first, A = Ao and all modes with k < ko (see equation 
(f80)l ) will go unstable, with the fastest-growing one given by 
kp ~ fco/\/2- Eventually the amplitude in this mode reaches 
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Figure 3. Left panel: spectrum of the magnetic fluctuation energy at three specific times {tun = 4, 16) during the evolution shown in 
figure [21 red short-dashed lines show the slope (equation l|96|l ): blue long-dashed lines show the firehose growth rate (see equation 
i79t ) for the instantaneous values of A -|- 2//3i at the times the spectra are plotted. Right panel: magnetic fluctuations in real space 
{Bix/Bq vs. z) at the same times as the spectra in the left panel. 



the level at which the back-reaction becomes important: ap- 
proximating equation (|87|) by 

1 d 



A(t) ~Ao + 



2v, dt 



(88) 



we find that the nonlinear contribution is comparable to 
|A() +2//3ij when (cf. equation ((30)l ): 



El 



7n 



2p, 
A. 



mfp 



(89) 



where 7max is the imaginary part of uip given by equation 
(|82p . Equation (|89p only gives a good estimate of the criti- 
cal amplitude if 3j^, is larger or not too much smaller than 
7max (collisions are sufficiently strong). If i/, ^ 7max (as a 
subsidiary limit within our e ordering), then a better ap- 
proximation than equation (|88|l is to replace the colhsional 
relaxation exponent in equation (|87p by unity, which gives 



E 



A, 



2 A 2 

3 + ^ 



(90) 



Once the nonlinear feedback becomes active, exponen- 
tial growth must cease and secular growth starts because the 
anisotropy must be kept close to marginal; using equation 
we find, to dominant order. 



2 



Ei^^ 



Ao + 



(91) 



This is valid regardless of which of the two estimates 
or (|90|l of the amplitude at the onset of nonlinearity was 
appropriate. This is because the effective growth rate asso- 
ciated with the secular growth decreases with time and so we 
will always eventually end up in the regime where the col- 
hsional relaxation exponent in equation l|87|l is faster than 
the magnetic energy growth and equation (|88p gives a good 
approximation of equation (|87|l . 

The evolution of the fiuctuation spectrum must be con- 
sistent with equation (|91|l . As the magnitude of the to- 
tal pressure anisotropy A approaches the marginal value, 
both the cutoff wavenumber fco(t) and the most unstable 
wavenumber fcp(t) decrease, as they can still be estimated 



by equations (|80|l and (|81|l with A = A(f). The modes whose 
growth has been thus switched off become oscillatory: from 
equation (|86[l . it is obvious that for k 3> ko{t), 



+ C2e 



(92) 



where ci and C2 are integration constants and c1 = C2 be- 
cause Al. — A-k (note that this oscillation of the ampli- 
tude is superimposed on the oscillation with the same fre- 
quency that was factored out in equations (|83H84p ). Since 
these modes oscillate in time at a rate that is much larger 
than the rate of change of the anisotropy, they no longer 
contribute to the feedback term in equation (|87|l . 

Thus, as the range of growing modes, peaked at kp{t) 
and cut off at ko{t), sweeps from large to small wavenumbers, 
they leave behind a spectrum of effectively passive oscilla- 
tions, whose amplitude no longer changes. Since there is no 
fixed special scale in the problem (except initial most unsta- 
ble wavenumber), one expects the evolution to be self-similar 
and the spectrum a power law. It is not hard to determine 
its exponent. Let \Ak\'^ ~ fc^". Since the total energy must 
grow linearly (equation (|9ip ) 



El 



t 



t' 



-l/{o 



(93) 



(this is valid if a > 1; the extra power of k comes from 
the integration over wavenumbers). On the other hand, for 
the fastest-growing mode, we must have, assuming secular 
growth, 

1 dA. 



Ak 



dt 



' 7n 



A + 



(94) 



where the last relation follows from equation (|82p . This gives 
us a prediction for the time evolution of the residual pres- 
sure anisotropy and, via equation ((81]), of the most unstable 
wavenumber (the infrared cutoff of the spectrum): 



-I 



1_ 

7t 



(95) 



The only way to reconcile equations (|93p and (|95p is to set 
Q = 3. Thus, we expect the one-dimensional firehose turbu- 
lence spectrum to scale as 
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(96) 



The secular growth of the firehose fluctuations will con- 
tinue until our asymptotic expansion becomes invalid, i.e., 
when the fluctuation amplitude is no longer 

smallEl 

From 

equation (|9ip . this happens at t ~ (uhIAq + 2/ I3i\)^^ ~ 
|7o|~^, where dimensions have been restored. This is the 
time scale of the large-scale dynamics. Thus, as we have al- 
ready explained in section 12.41 there is no saturation of the 
firehose fluctuations on any faster time scale. Unsurprisingly, 
at the same time as the fluctuation amplitude becomes large 
enough to break our ordering, the scale of the fluctuations 
also breaks the ordering: substituting the above time scale 
into equation ^5^, k\^\pi ~ (|7o|/r2i)^''^ ~ <^^i or fc||^mfp ~ 1, 
while our original ordering assumption was kupi ~ e, or 



fciiA 



mfp 



1/e (see equation (|54|l '). 



4-3.3 Numerical solution 

The firehose turbulence equation (|86|) is one-dimensional, so 
it is very easy to solve numerically; equation (|87|l is most 
conveniently solved in a differential form: 



(97) 



with the initial condition A(0) — Aq. Here we describe the 
results obtained from such a numerical calculation with the 
following parameters: 



Ao = -0.02, 



= O.OOOi. 



mfp 



(98) 



This means that the maximum wavenumber at which fire- 
hose fluctuations can be excited is fco — 0.28 (equation (|80p : 
see figure We solve equation (|86|) for 1024 wavenum- 
bers in a periodic domain of size Amfp, so the smallest 
and the largest wavenumbers are (still normalized to pi) 
fcniin = 27rpi/Amfp — 0.00063 and fcmax — 0.32. The initial 
conditions are random amplitudes in each wavenumber (sat- 
isfying the reality condition A^k = AX)- Note that with the 
parameters (|98[) . our ordering parameter is e ~ 0.1, so we 
have chosen a spatial scale separation between collisions and 
the Larmor motion that substantially exceeds 1/e^ formally 
mandated by our ordering (section l3.4p . This does not break 
anything and is in fact more realistic for the physical param- 
eters in weakly coUisional plasmas of interest (section l2.ip . ft 
also widens the scale interval available to the firehose turbu- 
lence spectrum and ensures that even deep in the nonlinear 
regime, when the wavenumber of the firehose fluctuations 



Note that while the amplitude grows and thus eventually 
breaks the ordering introduced in scction[3]4] the stability param- 
eter |A + 2//3i| decreases, so the approximation of small Larmor 
radius gets quantitatively better with the growth of the firehose 
fluctuations moving to larger scales (equation (|8TJ) — equiva- 
lently, our ordering of pi introduced in section FS. 41 ieg nation (|56j) 
is quantitatively better satisfied. In fact, we could have chosen 
to construct our entire asymptotic theory by expanding close to 
marginal stability and so ordering everything with respect to the 
small parameter defined as e = |A -|- 2/(9; [^/^ instead of equation 
1451 1 (this is t he route followed in an analogous mirro r instability 
calculation bv lRincon. Schekochihin fc Cowievll2010h . 



drops substantially, there is still a healthy scale separation 
between them and the coUisional dynamics. 

The evolution of the total magnetic energy, \Bi\^/Bq = 
X]fe l^fePi is shown in figure [2] (left panel). Initially it grows 
exponentially at the (normalized) rate 7max = Imcjp (see 
equation (|79|l : this part of the evolution is trivial and so 
not shown). The exponential growth is followed by a secu- 
lar, linear in time, growth of the energy in accordance with 
equation (|9ip . The energy at which this nonlinear regime 
starts is closer to the estimate given by equation (|90|l than 
by equation (|89p because, as discussed above, we have taken 
a very small value of Note that in this and all subse- 
quent figures, we have normalized time using the collision 
frequency vu, not the cyclotron frequency Qi — this is indi- 
cated explicitly in the figures and should cause no confusion 
to an attentive reader. 

Figure [5] (right panel) shows the time evolution of the 
instability parameter A -|- 2/ Pi. As expected, it is tending 
to the marginal stability value (zero). The inset shows that 
this approach to zero is consistent with the 1/t prediction 
(equation (|95j)E3 

The evolution of the spectrum of firehose fluctuations is 
illustrated by flgure |3] (left panel) . As anticipated in section 
14.3.21 the spectral peak moves to smaller wavenumbers in the 
nonlinear regime. The spectrum extending from this moving 
peak to the original wavenumber of the fastest linear growth 
(fcp = 0.2; see equation (|8ip ) is statistically stationary and 
consistent with the fc~^ power law predicted by equation 
(|96p . The instantaneous flrehose growth rate is overplotted 
on the spectra in figure |3] (left panel) and confirms that the 
position of the spectral peak closely follows the wavenumber 
of the fastest instantaneous growth of the firehose instability. 

Figure [3] (right panel) shows snapshots of one of the 
components {Bix) of the perturbed magnetic field corre- 
sponding to the spectra in figure [3] (left panel). The emer- 
gence of increasingly larger-scale fluctuations is manifest. 
Perhaps a better illustration of this real-space evolution of 
the firehose turbulence is figure |4] (left panel), which is the 
space-time contour plot for the middle fifth of the domain. 



4.4 Implications for momentum transport 

Substituting the second-order pressure tensor calculated in 
section [T71 into the large-scale momentum equation (|7ip . we 
get 



niinoi- 



dt 



-Vp + V ■ 



Poiboba 



-I 



(99) 



where, in the absence of density and temperature gradients, 
the total isotropic pressure p = P2i + /St!" is set by the con- 
dition V-iioi — Or^l while the pressure anisotropy A is given 
by equation H73p . A remarkable feature of equation (|99p is 
that all of the effects of the magnetic field appear in the term 



The oscillations seen in the figure are no t a numerical arte- 
fact. T hey are due to oscillatory transients — ISchekochihin et al.l 
1I2OO8I ') derived those analytically for a solution with only one 
Fourier mode. 

More generally, p adjusts in such a way as to reconcile the pres- 
sure balance with the continuity and heat conduction equations; 
see Appendix |A2. 121 
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Figure 4. Left panel: contour plot of the time evolution of the firehose turbulence for the numerical solution discussed in section 14.3.31 
— the vertical axis is space (the middle fifth of our entire periodic domain), horizontal axis is time, the colours respresent the value 
of Bix/Bq. Right panel: the same plot for the gyrothermal turbulence discussed in section [5.3.31 Note that, as explained in the text, 
the firehose turbulence exhibits a gradual coarsening of the dominant structure with time, while the gyrothermal turbulence ends up 
dominated by a single scale. 



proportional to A + 2//3i, which is precisely the instability 
parameter that the small-scale firehose turbulence described 
in section f43l contrives to make vanish. In the marginal state 
that results, the tension force (the 2 /Pi term) is almost en- 
tirely cancelled by the combined pressure anisotropy due to 
large- and small-scale fields. This suggests that in regions of 
the plasma where the firehose is triggered (i.e., where the 
magnetic field is locally decreased by the plasma motion), 
the plasma motions become effectively hydrodynamic, with 
magnetic-field lines unable to resist bending by the flows. 



Since the cancellation of the second term in equation 
(|99p by the firehose turbulence also effectively removes the 
(parallel) viscosity of the plasma, these hydrodynamic mo- 
tions are not dissipated. In a turbulent situation, this should 
enable a cascade to ever smaller scales. Obviously, once this 
happens, the original motion that caused the negative pres- 
sure anisotropy to develop is supplanted by other, faster 
motions on smaller scales. The theory developed above even- 
tually breaks down because the scale separation that formed 
the basis of our asymptotic expansion is compromised: while 
the fiuid motions penetrate to smaller scales, the firehose 
fluctuations move to larger scales (see section fOjl . 



Note also that the fluid motions produced by the tur- 
bulent cascade can give rise to both positive and negative 
pressure anisotropies — and so, to have a full description 
of their further evolution, we must know the effect on mo- 
mentum transport not just of the firehose but also of the 
mirror and other instabilities triggered by positive pres- 
sure anisotropies (locally increasing magnetic field strength). 
Th is is still work in progress (the mirror case is considered 
by iRincon. Schekochihin fc CowlevI I2OI0I '). Another impor- 
tant adjustment to the viscous-stress reduction argument 
above has to do with the modification of the firehose insta- 
bility by the parallel ion heat fluxes — we now proceed to 
investigate this. 



5 GYROTHERMAL TURBULENCE 

5.1 Firehose turbulence equation with heat fluxes 

As we briefly mentioned in section 13.71 allowing a non- 
zero ion temperature gradient along the unperturbed mag- 
netic fleld leads to substantial modiflcations. These are of 
two kinds. First, as shown in Appendix IA2.13I the pres- 
sure anisotropy Aq caused by the large-scale dynamics con- 
tains contributions from the coUisional parallel heat fluxes 
(proportional to 60 ■ VToi) and from compressive motions 
(as we pointed out in footnote (9] the presence of a tem- 
perature gradient automatically implies a density gradient 
as well because of the requirement that pressure balance 
should be maintained; see equation HA36|) and Appendix 
IA2.12|I . Instead of equation ((74]), valid in the incompressible 
case, we must use the more general equation (|A63|l . This, 
however, does not change much: the unstable flrehose fluc- 
tuations will grow in the manner described in section 14.31 
first exponentially, then secularly, to compensate whatever 
pressure anisotropy is set up by the large-scale dynamics. 
The only change is the physical interpretation of the origin 
of the pressure anisotropy: as long as ion temperature gra- 
dients are present, the anisotropy is not tied exclusively to 
the change in the magnetic field. Physically, the heat-flux 
contributions to the anisotropy have to do with the fact 
that "parallel" and "perpendicular" heat flows along the 
magnetic-fleld lines somewhat differently and so imbalances 
between p± and py can occur — this can be seen already 
from the CGL equations (see Appendix IA2.15|) . 

The second heat-flux-related modiflcation of the theory 
developed thus far is more serious. It involves an additional 
contribution to the FLR term in the third-order pressure 
tensor (equation (|75|) ') and, therefore, to the flrehose turbu- 
lence equation (|77p . This contribution was derived in Ap- 
pendix |X2T4l but suppressed in our previous discussion. It 
is given by equation (|A64|) and consequently equation (|77|) 
now reads 

V||[(^A + -JB, 
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Figure 5. Left panel: frequencies (thin lines) and growth rates (bold lines) of the unstable firehose modes (red/solid: the "+" mode; 
blue/dashed: the "— " mode) given by equation 1)103^ : the parameters here are A + 2//9i = —0.01 and Ft = 0.02, so the instability param- 
eter is A = 0.0054 (equation I I104I I'). Right panel: same, but for A + 2/l3i = 0.00075, so A = 0.000025 (close to marginal stability); dotted 
vertical lines indicate the wavcnumber of fastest growth kp = 0.085 (equation H109p ) and the dotted horizontal lines the corresponding 
maximum growth rate 'ymax — 

Imwp = 0.0004 (equation i fTTOll '). 
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We have introduced a dimensionless parameter measuring 
the magnitude of the parallel heat fluxP^ 



1 vthi bo ■ VToi 1 A 
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(101) 



where It is the parallel length scale of the ion temperature 
variation. We see that the functional form of the firehose 
turbulence equation is changed. We now proceed to study 
the effect of this change. 



5.2 Linear theory: the gyrothermal instability 

The linear dispersion relation for equation (|100p is 



.(102) 



Like in the case of equation (|78|l . there are four roots of 
which two are potentially unstable: 



(103) 



We stress that we are discussing the efl'ect of the ion heat 
flux as the electrons are assumed isothermal at the scales we are 
considering (see Appendix lAljl . We also stress that these heat- 
flux effects enter through the FLR terms in th e plasma p r essure 
tensor and are absent in, e.g., the lowest-order iBraeinskiil lll965l ) 
equations. 



where k = k^^pi. Instability occurs at wavenumbers for which 
the expression under the square root is positive. There is an 
interval of such unstable wavenumbers if and only if 

A-r|-i(A + l)>0. (104) 

If this condition is satisfied, the "-|-" mode is unstable for 

- 4 (Lt + Va) < fc < -4 (Lt - Va) , (105) 



(106) 



and the "— " mode for 

4 (Tt - VlCj < fc < 4 (^Ft + Va) . 

where we have assumed, without loss of generality, that 
Ft > 0. When A + 2//3i < 0, these two intervals intersect, so 

all modes with |fc| < ko ~ are unstable (oth- 

ers are pure propagating waves). When A + 2//3i > 0, the 
intervals are separated and there is an interval of stability 
at long wavelengths, viz.. 

What is remarkable about all this is that not only the 
stability conditions and specific expressions for the firehose 
growth rate are modified by heat flux, but the presence of 
the heat flux allows for instability even when firehose is sta- 
ble, A + 2/ Pi > (but positive pressure anisotropy not 
too large and /3i not too small, subject to equation (|104|) ). 
This instability, called the gyrothermal instability ( GTI), 
leads to the growth of Alfvenically polarized fluctuations 
in the parameter regime in which they are otherwise stable 
ISchekochihin et al.ll201ol )F1 



^"^ Note that for A — l/ft > 0, the mirror mode is unsta- 
ble as well, but it involves growth of compressive fluctuations, 
SBn 3> <5iJxi at highly transverse wavenumbers fcii <C fcx (see, 
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Figure 6. Left panel: evolution of the magnetie energy |Bj^p/Bg = l^fcP with time in a numerical solution of equations l|112|l and 
II113II with parameters I I118I I: the red dotted line shows the nonlinear asymptotic given by equation II114II : this figure is the GTI analog 
of figure [2] (left panel) . Right panel: evolution of the instability parameter A (pink) and the pressure anisotropy parameter A + 2//?^ 
(black) in the same numerical solution. Inset: log-log plot of the evolution of |A|; the black line shows the slope corresponding to 
(see equation l|116|l ). 



The formulae for the wavenumber of the fastest-growing 
mode and the maximum growth rate for the combined 
firehose-GTI are straightforward to write down. As always 
with such formulae, they are not particularly illuminating in 
the general case, but are interesting in various asymptotic 
limits. When the firehose instability parameter A-|-2//3i < 
and its magnitude is much larger than Fy, the effect of 
the heat flux is a small correction to the firehose insta- 
bility already described in section 14.21 Conversely, when 
|A -f 2//3i| < r^, the GTI is dominant and, for the fastest 
growing mode, 

fcp ~ =F6rT, 

cjp ~ 9r|. ( ±1 + 



(107) 
(108) 



where Up is normalized to Qi . Finally, close to the marginal 
state, A — -f 0, we have 



fcp ~ =F4rT ( 1 + p2 I ' 



cjp ~ 4r^ ±1 + 



A 

T 

^^/A 



(109) 
(110) 



Note that, unlike the firehose, the GTI has a definite pre- 
ferred wavenumber that does not change as marginal stabil- 
ity is approached. 

Figure [S] (left panel) shows the dependence of the fre- 
quencies and growth rates of the two unstable modes on 
wavenumber for a set of parameters for which the instabil- 
ity is a hybrid of firehose and GTI (these are the parame- 
ters used in the numerical solution of section [5. 3. 3|) . Figure 
[S] (right panel) shows the same for the case in which the 
firehose is stable (A -I- 2//3i > 0) and that is very close to 
marginal stability: we see that the instability only exists in 



e.g.. lHellingei1l2007l 'l. while Alfvenic fluctuations are not affected 
by it to lowest order in the instabi lity parameter A — l//3i 
llRincon, Schekochihin fc Cowlevll201Cll ). 



the immediate neighbourhood of the last unstable wavenum- 
ber given by equation (|109[l . 



5.3 Nonlinear evolution and spectrum 

5.3.1 Firehose-GTI turbulence equation in scalar form 



As happened in section 14.3.11 equation (|100|l can be re- 
duced to one equation for a scalar field, although it is now a 
slightly more complicated transformation. Let us again non- 
dimensionalize time and space according to equation 
and introduce new fields A'^{t) as follows: 
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With the ansatz (|llip . equation (|100|l becomes 
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A{t) = Ao + 
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It is now manifest how the dispersion relation (|103[l emerges 
from equation (|112|l . Unlike in the case of pure firehose tur- 
bulence (Ft = 0), the evolution of the mode now depends 
on the sign of its real frequency — that is why we have two 
scalar equations. However, these equations have a symme- 
try: if we arrange initially that A'^ = AZj. (which we can 
always do by an appropriate choice of the phases (^fe), then 
this relation will continue to be satisfied at later times. This 
also means that A^ are real because, in order for Bi" to be a 
real field, we must have (from equation (jllip ) (A'^)* = AZ^. 
(we assume the phases satisfy (j>k = (t>-k). The conclusion is 
that it is enough to solve just one of the two equations (|112|l 
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Figure 7. Left panel: spectrum of the magnetic fluctuation energy at three specific times {tun = lj4, 16) during the evolution shown in 
figure [6l red short-dashed lines show the slope; blue long-dashed lines show the firehose/GTI growth rate (see section f5.2ll for the 
instantaneous values of A + 2//3i at the times the spectra are plotted. Right panel: magnetic fluctuations in real space {Bi^/Bo vs. z) 
at the same times as the spectra in the left panel. This figure is the GTI analog of figure |3] 



— either for the -I- or the — mode. The total energies of the 
two modes that enter equation (|113|l are equal0 



5.3.2 Qualitative picture 

The evolution of the firehose-GTI turbulence is easy to 
predict arguing along the same lines as we did in section 
14.3.21 Let us consider the case when initially the pressure 
anisotropy is negative and — (Ao -I- "2/ Pi) S> 2r|., i.e., the 
instability parameter Ao > (given by equation (|104|) with 
A = Ao). In this regime, the heat flux does not matter and 
the evolution proceeds as in the case of the firehose turbu- 
lence: magnetic fluctuations grow and eventually the nonlin- 
ear feedback in equation (|113p starts giving an appreciable 
positive contribution to the pressure anisotropy (estimates 
(|89p and (|90|) for the fluctuation amplitude at which this 
happens are still valid). A k^^ spectrum will then form, with 
the infrared cutoff (wavenumber of maximum growth) mov- 
ing to larger scales and |A-|-2//3i| decreasing (i.e., A-f 2//3i 
increasing and thus becoming less negative). 

The evolution of the gyrothermal fluctuations starts to 
differ from the pure firehose case after |A -I- 2//3i| becomes 
comparable to V%. The GTI is now the dominant insta- 
bifity mechanism. Since the fiuctuations continue growing, 
A -|- 2/ Pi continues to increase and will become positive, 
tending eventually to 2V%, so as to push the instability 
parameter A (equation (|104|) ') to zero and the GTI to its 
marginal state. As A — > +0, the growth is concentrated in 
a shrinking neighbourhood of the wavenumber ftp = 4rT 
(see equation (|109p ). This means that the spectrum stops 
spreading towards lower wavenumbers and its infrared cut- 
off stabilizes at fcp. All the growth of magnetic energy is now 
provided by the growth of the one mode associated with fcp, 
which will soon tower over the rest of the spectrum. 



The same approach could have been taken in section r4.3.1l in- 
stead of solving equation I I112I1 for a complex function Aj. subject 
to = A_j., we could have solved for one of two real functions 
subject to = A"^.. The magnetic field is then recovered 
via equation milll . 



The growth is still secular: using equation H113P and 
the marginaUty condition A = 0, we find to dominant order, 
analogously to equation H91|) . 

^ > ^lAfcl^ ~4Ao;/.i. (114) 
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Finally, we can calculate the evolution of the residual A. 
Analogously to equation (|94p . the growing mode satisfies 

1 dA. 1 



dt 

where we used equation (|110|) for 7n 

A' 



Therefore, 
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(115) 



(116) 



As in the case of the firehose turbulence, the secular 
growth will continue until the fiuctuation amplitude is no 
longer small: t ~ (i/uAo)'^ ~ l7o|~^ (time scale of the large- 
scale dynamics). The key difference from the pure firehose 
case is that the fiuctuations are now stuck at a microscopic 
spatial scale given by equation (|109|l : restoring dimensions 
and using equation (|10ip . the corresponding wavenumber is 

Amfp 



kupi 



h 



(117) 



(this scale is collisionless, fcyAmfp ^ 1, provided It <^ 
-^mfp/pii for galaxy clusters, this is always true as is easy 
to ascertain by using the numbers from section ETTj) . Thus, 
the gyrothermal turbulence is essentially one-scale, in the 
sense that fiuctuations at this one scale become energeti- 
cally dominant as marginal stability is approached at late 
stages of the nonlinear evolution. 

5.3.3 Numerical solution 

We have solved equations (|112p and (|113p in a manner com- 
pletely analogous to that described in section [4.3.31 The pa- 
rameters we used are 



Ao = -0.02, 
pi 



u. 



X 



P^ 
0.0001 



0.01, Ft = 0.02, 



mfp 



(118) 
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This implies that the instability parameter in the linear 
regime is A = 0.0054 (equation (|104|) ) and so the maximum 
unstable wavenumber is ko — 4(rT + \/A) — 0.37 (equations 
HlOSp and (|106|) : see figure [5] (left panel)). Our numerical 
solution now has 2048 wavenumbers, so kmin — 0.00063 and 

fcmax — 0.64. 

As expected, the evolution of the total magnetic energy 
is similar to the case of pure firehose turbulence discussed in 
section [4.3.3l exponential, then secular growth (see equation 
H114P ') — this is shown in figure[S](left panel). The evolution 
of the instability parameter A (equation (|104|l ') towards its 
zero marginal value is given in figure [S] (right panel) . The 
inset shows that this approach to zero is consistent with 
the 1/t^ prediction (equation (|116|) ). Also shown in figure 
[H] (right panel) is the evolution of the pressure anisotropy 
parameter A + 2/ Pi, which for the pure firehose used to be 
the instability parameter. Since A — ^ 0, it should tend to 
2rl = 0.0008 and it indeed does. 

Finally, figure [7] (left panel) illustrates the evolution of 
the spectrum of firehose/gyrothermal fluctuations. It fol- 
lows the scenario outlined in section 15.3.21 At first it is 
similar to the firehose turbulence spectrum with the spec- 
tral peak moving towards larger scales leaving behind a 

spectrum. As the wavenumber of fastest growth kp 
approaches the value corresponding to the near-marginal 
GTI, kp = 4rT = 0.08 (see equation (|109l) ). the peak stays 
there and continues growing, eventually dominating all other 
modes. The emergence of a one-scale sea of gyrothermal fluc- 
tuations is further illustrated by figure |4] (right panel), which 
shows what these fluctuations look like in real space as time 
progresses. The difference between them and the pure flre- 
hose fluctuations in figure |4] (left panel) is manifest: the gy- 
rothermal ones stay at the same scale while the firehose ones 
become larger-scale as time progresses. 



5.4 Implications for momentum and heat 
transport 

Let us now revisit the discussion of the effect of plasma 
instabilities on the momentum transport modification at- 
tempted for the pure firehose in section [J?] As before, the 
combined large-scale viscous and Maxwell stress is contained 
in the second term on the right-hand side of equation (|99p . 
However, with parallel ion heat fiuxes present, the nonlin- 
ear evolution of the GTI pushes the quantity A -f 2 /Pi not 
to zero but to a positive value 2T'^, corresponding to the 
marginal state A = (equation (|114|) '). Since any smaller 
value of A + 2 /Pi is GTI unstable, this leads to a curious 
conclusion that the momentum transport is now effectively 
determined by the ion heat fiux: 

niinoi^^^^ = -Vp+ V- (poibobo2V^ 



dt 



b() bo- 



no; ( bo • VTo, 



rriiul 



(119) 



where we used equation (|101|l for Ft. This equation has to 
be supplemented with the transport and pressure-balance 
equations for no; , To; and p as explained in Appendix lA2. 12l 
Equation (|119p probably merits a careful study (which 
is outside the scope of this paper), but we would like to 



accompany it with a very important caveat. Since pressure 
anisotropy in the nonlinear state of the GTI can be positive, 
other plasma instabilities may be triggered. Thus, if A > 
1/ft i.e.,ffF;^ > 3/ (2ft), the pi asma will be mirror unstable 
(see IHellingerl |2007| . and references therein) . The magnetic 
fluctuations that the mirror instability produces are different 
from the GTI both in polarization (5-B|| , not SB±) and scale 
(fc_L > fe||, k_LP^ ~ (A - l//3.)'/^ ~ A 1/p, for the 
mirror, whereas for the GTI we had k± — 0, ~ 4Ft). 
How they saturate and what they do to the effective pres- 
sure anisotropy is a matter under active current investiga- 
tion ICalifano ct al. 2008: Istomin. Pokhotelov fc BalikhinI 
l2009l : iRincon. Schekochihin fc CowlevI bOloD — and it is 



completely unknown how mirror and gyrothermal fluctua- 
tions might coexist. 

The key question is whether the pressure anisotropy 
will be set by the GTI or the mirror marginal condition 
and if it is set by the latter (A = 1/Pi, as, e.g., seems 
to be indicated by the solar wind data; see iHellinger et al.l 
l2006l : iBale et al.1 [20091 ). then whether a turbulent plasma 
has a way of suppressing the GTI by adjusting not the pres- 
sure anisotropy, but the heat flux to the marginal condi- 
tion: r'^ = 3/(2^0- This raises the possibility that not only 
the pressure anisotropy but also the (ion) heat fluxes are 
determined by the marginal stability conditions of the fire- 
hose/GTI and mirror. Thus, plasma instabilities may be the 
crucial factor in setting both the momentum and heat trans- 
port properties of a weakly collisional plasma. We stress, 
however, that under the assumptions adopted in this paper, 
we have not produced a nonlinear mechanism for changing 
the ion heat flux and this remains a subject for future work. 



6 DISCUSSION AND COMPARISON WITH 
PREVIOUS WORK 

6.1 Marginal stability via particle scattering or 
via changing field structure? 

It is not in itself particularly surprising that the nonlin- 
ear effect of an instability driven by pressure anisotropy 
is to produce fluctuations that effectively pin this pres- 
sure anisotropy at a value corresponding to marginal sta- 
bility. Bes ides having direct observational support in the so- 
lar wind (iGarv et al.l 200ll: Kasper. Lazarus fc Gary! 20021: 



Marsch. Ao fc Tull20q4l : lHeUinger et al.ll2006l : lMatteini et all 



20071 : iBale et al.ll2009^ ~ it makes sense as a fundamental the- 
oretical expectation (Le Chatelier's principle). One may be 
tempted to proceed to another, seemingly as reasonable, the- 
oretical expectation that the mechanism for achieving this 
marginal state must be pitch-angle scattering of particles by 
the fluctuations leading to isotropization of pressure. While 
indeed physically reasonable, this is, however, not an in- 
evitable conclusion. As we have shown above, particle scat- 
tering is, in fact, not the way the k± — firehose fluctuations 
make pressure anisotropy marginal (under the ordering as- 
sumptions we have adopted). Instead, the marginal state is 
achieved via a modification of the structure of the magnetic 
field: namely, secular growth of the microscale fluctuations 
cancels on average the decrease in the mean field that pro- 
duced the pressure anisotropy thus pushing the latter to its 
marginal value. This was explained on an intuitive level in 
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section 12.41 and the subsequent analytically rigorous devel- 
opments showed that intuition to be correct. 

Considering this result, we must recognize it as physi- 
cally reasonable on the following grounds. A particle trav- 
elling in a magnetic field will traverse a fluctuation with a 
given fc|| over time ~ l/fc||Dthi. This time is much longer than 
the ion cyclotron period if fcyUthi ^ ^i, or, equivalanetly, if 
fc||Pi <^ 1. If this condition is satisfied and if the frequency 
of the fluctuation uj <^ Qi, the fluctuation cannot change 
the first adiabatic invariant fi — v\/2B of the particle, so 
there cannot be very much pitch-angle scattering. In our 
calculation, as the pressure anisotropy (or, more precisely, 
the instability parameters A + 2//3i and Ft) were small, the 
parallel scale of the fluctuations generated by the k± = Q 
firehose or gyrothermal instabilities was substantially larger 
than the Larmor scale (see sections 14.21 and I5.2[l and, in 
the case of the firehose, it increased further in the nonlin- 
ear regime (see section l4.3.2p . Thus, k\\pi <C 1 was satisfied 
at all times (as was uj <^ Qi), the plasma remained magne- 
tized and pitch-angle scattering ineffective, so the rearrange- 
ment of the field structure was the only device available to 
the system to counteract the pressure anisotropy drive. It 
is possible that the oblique firehose (which is much harder 
to treat analytically than the parallel one) might produce 
fiuctuations at the ion Larmor scale, so particle scattering 
by firehose fiuctuations is not completely ruled out, but it 
certainly does not happen for the fc^ = case to which we 
have limited the scope of the present investigation. 

How important is it to know whether particle scattering 
is present? Re c ently, in the context of accretion-disc physics, 
[Sh arma et all (|2006l . 120071 ') proposed an ad hoc closure for 
numerical simulations, constraining the pressure anisotropy 
to lie within the marginal stability boundaries via artificial 
dissipation in the pressure equations (the CGL equations 
given in Appendix I A2. 15| l . They argued that this was justi- 
fied if it could be shown microphysically that plasma insta- 
bilities (in their case, ion cyclotron and firehose) produced 
fiuctuations at the ion Larmor scale, where pitch-angle scat- 
tering of particles off the fiuctuation "foam" isotropized 
pressurel^ As we have explained, our results for the par- 
allel {k± = 0) firehose do not support this picture. How- 
ever, it is not obvious that the validity of a closure based 
on the average pressure anisotropy being maintained at the 
marginal level must be predicated on the presence of parti- 
cle scattering. As we have shown above, a sea of secularly 
growing magnetic fiuctuations far above the Larmor scale 
can produce the same effect. This, of course, does not ex- 
cuse us from having to find the right microphysical theory 
for pressure isotropization if we are ever to have anything 
more than a plausible closure imposed by fiat. 

One example of a context in which the presence or 
absence of scattering matters greatly is viscous heating of 
the plasma. The heating depends both on the pressure 
ani sotropy and on t he collision frequency (it is oc unA^; 
see iKunz et al.ll201ll ). so, in order to calculate it correctly, 
we must know whether only the pressure anisotropy or also 



Tlie same view was taken bv lSchekochihin fc Cowlevi | |2006| ) in 
tlieir m odel of the dynam o action in a weakly collisional plasma 
and by iBale et al.l ||2009| ) in interpreting their measurements of 
marginal pressure anisotropics in the solar wind. 



the (effective) collision frequency is modified by the fire- 
hose fluctu ations. Assuming t he Coulomb collision frequency 
unchanged. IKuuz et al.ir201ll recently proposed a thermally 
stable heating mechanism for galaxy clusters (see a further 
short discussion in section [7.3.ip . If microphysically justified 
in the most general case (i.e., not only for the parallel fire- 
hose, but also the oblique one, the mirror instability, etc.), 
this represents significant progress. Thus, having a detailed 
microphysical theory does make a difference not only for 
the analytical strength of the subject but also for explaining 
astronomically observed realities. 



6.2 Quasilinear theories 

The theory developed above is basically quasilinear in that 
the fiuctuation amplitude is assumed small and it is found 
that such small fiuctuations can drive the instability to a 
marginal state. There have been a number of quasilinear 
treat m ents of the firehose instabi lity fehapiro fc Shcvchcnkd 
19641: iKennel fc SagdeCTi 1 19671: iDavidson fc Volk ,196^ : 



Gary fc Feldmanlll978l : iQuest fc Shapirdll996l ). so it is per- 
haps useful to explain why they do not obtain similar results. 

The approach in such theories is to consider a collision- 
less plasma with some initial distribution that has a negative 
pressure anisotropy (let us call it Ao < 0) and work out how 
it relaxes. The result is that a fiuctuation level builds up. 



witlO 



(120) 



which is small when the instability parameter Aq -f 2//3i is 
small. This saturated fiuctuation level suffices to marginalize 
the instability. This result is easily recovered in our theory 
if we formally set Ua = in equation (|73p . This gives 



A{t) = Ao + 



3\BUt)l 
2 



(121) 



and assuming saturation in the marginal state A = — 2//?^, 
we recover e quation (I120II. The classic wo rk where this was 
first done is IShapiro fc Shevchenkd (|l964l 'l (a more detailed 
comparison is provided in Appendix IA2. 16) ) . We stress that 
in their calculation the saturation of the pressure anisotropy 
at the marginal level is not due to particle scattering any 
more than it was in ours because the fiuctuations still have 
<^ 1, <jj <C 0.i and so conserve (see section 16. ip . 
The internal energy stored in the pressure anisotropy is 
transferred into magnetic fiuctuations until the pressure 
anisotropy is marginal. The magnetic fiuctuations then per- 
sist because under the adopted appriximations there is no 
dissipation of the magnetic field. 

The difference in our approach is to include weak col- 
lisions and consider the case when the pressure anisotropy 
is constantly driven by the large-scale dynamics (which is 
physically where it comes from; see section [2^ . The steady 



20 iHalll iligsiT) argues qualitatively for a similar saturation level, 
but due to trapping of particles in firehose fluctuations. As the 
systematic kinetic calculation presented above shows, the trap- 
ping eff'ect does not play a role at these amplitudes, at least not 
under the assumptions we adopted {k^_ = and relatively high 
collision frequency). 
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level of the anisotropy is then set by the competition be- 
tween collisions and the drive (equation (|74l) ') and to off- 
set this anisotropy and keep the instability marginal, the 
fluctuation level has to keep growing secularly rather than 
stay constant (the earlier quasilinear theories can then be 
interpreted to describe correctly what happens before one 
collision time has elapsed). 



6.3 Driven anisotropy in a collisionless plasma 

It is interesting to inquire what would happen if the 
anisotropy were driven (rather than just initially imposed) 
but collisions not strong enough to balance the drive and 
impose a steady anisotropy. Formally speaking, our the- 
ory breaks down in this case because the equilibrium dis- 
tribution cannot be proved Maxwellian. However, that is a 
technical issue and one could, in fact, reformulate our the- 
ory as a near-marginal expansion in the instability parame- 
ter |A -I- 2//3i|. We expect that equation (|A60[l (or, equiva- 
lently, equation (|A67|l '). with the coUisional relaxation term 
removed, would still describe the evolution of the anisotropy: 



dt 2dt Bl ' 



(122) 



where 70 (assumed negative) is the drive — it contains all 
the terms in equation (|A60[l due the large-scale dynamics. 
Under these conditions, the driven part of anisotropy is con- 
stantly increasing and so again the fluctuations will have to 
grow secularly in order to keep it at the marginal level: 
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This will, of course, break down once the fluctuation level is 
no longer small or collisions catch up. 

6.4 Otlier nonlinear theories and simulations 

There exist a number of numerical studies of the nonlinear 
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adopted the same relaxation-of-initial-anisotropy approach 
as the quasilinear theories discussed above and the results 
they report are broadly consistent in that the magnetic 
fluctuation energy saturates at a level scaling with the size 
of the initial anisotropy (see e quation (I120II) . A no table 
exception is the recent work of IMatteini et al] (|2006l ) who 
consider the anisotropy driven by the expansion of the solar 
wind — the fluctuation levels they see are probably well 
described by equation (|123|) . 

The spectrum of the flrehose fluctuations has not pre- 
viously been addressed analytically, but, perhaps vaguely 
in agreement with the results of section 14.31 some of the 
numerical evidence does point to the growing predomi- 
nance of smaller wavenumbers in the nonli near regime — as 
the anisotropy approaches marginal level (jOuest fc Shapirq 
1 19961 : IMatteini etaLlbOOd ) . 

Finally, to our knowledge, the effect of heat fluxes on the 
nonlinear behaviour of the flrehose instability, studied in sec- 
tionO has not been specifically considered before. Note that 



althoug h the heat fluxes are presen t in the numerical simula- 
tions of ISharma et al.l (|2006t 120071 ), their momentum equa- 
tion does not have the gyroviscous and gyrothermal terms 
that regularize the firehose at small scales and give rise to 
the gyrothermal instability. The appropriate modifica t ion to 
the fiuid equations suggested bv ISchekochihin et al.l (|2010l ) 
should in principle enable one to study the spectrum of the 
firehose and GTI fiuctuations numerically. 



7 ASTROPHYSICAL IMPLICATIONS 
7.1 Solar wind 

Much of the observational evidence about the firehose insta- 
bility comes from the measurements of pressure anisotropics 
and fiuctuation levels in the solar windj^ Since the wind is 
expanding, both the local density and the local magnetic- 
field strength are dropping, so one expects a negative pres- 
sure anisotropy to develop: this can be described by equa- 
tion p22|l . where the drive is roughly 70 ~ —Vs^/R (so- 
lar wind speed divided by the distance from the Sun) 
and the coUisional relaxation is neg lectedEI The evi- 
dence for this trend, negative pressure anisotropy devel- 
oping with increasin g distance from the Sun, is given by 
IMatteini et al.l |2003); a number of other papers also docu- 
ment the fact that the measured pressure anisotropics are 
bound ed from below by the fir ehose marginal stability con- 
dition ([Kasper. Lazarus fc Garv„2002 : HcUingcr ct al.,i2006l : 
iBale et al.ll200g " 



3^ 

Bale et al.1 (|2009t ) found increased levels of ion-Larmor- 



scale magnetic fluctuations close to this stability bound- 
ary — presumably due to the flrehose instability. There 
are also indications of an injection of energy int o paral- 
lel wavenumbers just ab ove the ion Larmor scale jPodestal 
I2OO9I : IWicks et al.l I2OIOI ) — again, conceivably by the flre- 
hose instability. It is unclear how these firehose fiuctuations 
coexist with the solar wind inertial- and dissipation-range 
turbulence — as our theoretical understanding of this turbu- 
lence is still largel y based on assuming isotr opic equilibrium 
distributions (see ISchekochihin et al.l I2OO9I . and references 
therein) . This is one of the contexts in which the absence of 
a complete microphysical theory of the firehose turbulence 
and its effect on the plasma motions is particularly acutely 
feh. 



7.2 Accretion discs 

Another such astrophysical context is hot accretion 
flows, in which the theoretical modelling of the long- 
standing problem of the angular momentum trans- 



There are also some measurements indicating the pres- 
ence of firehose fluctuations in the Earth's m agnetotail (see 
iHorton. Xu fc Wond |2004 iHorton et al] |2004 and references 
therein). 

^■^ Note, however, that although the mean free path in the solar 
wind is roughly comparable to 1 AU, one does find in the solar 
wind a strong correlatio n between press ure anisotropy and the 
estimated coUisional age llBale et al.l2009h . so modelling the solar- 
wind plasma as completely collisionless is possibly less valid than 
it might appear. 



Theory of firehose & gyrothermal instabilities 21 



port and radiative efficiency or inefficiency of the ac- 
cretion has taken a ne w turn with the introduction 
of p r essure anisotropics d Quatacrt. Borland fc ffammet; 
20021: ISharma. Hammett fc Quata cr t ,2003; Balbus ,200 



ion 

qI; 



Islani fc BalbusI l2005l: ISharnia et al.1 i2006l. |2007| ). The 
merical model of ISharma et al .1 (|2006l . 12007^ consisted of 
a closure that pinned down the pressure anisotropics at 
marginal stability via artificial dissipation terms. As we ex- 
plained in section 16.11 their assumption of microscale fiuc- 
tuations scattering particles is not bourne out by the theory 
developed above for the parallel {k± — 0) firehose fluctua- 
tions, although it is not excluded for other instabilities and, 
in any event, a closure based on marginal stability is prob- 
ably a sensible choice. 

A key remaining unknown here is the fate of the mi- 
croscale fluctuations over long (transport) time scales and 
the eventual structure of the tangled magnetic fleld that re- 
sults — a crucial question for accretion theories because they 
require knowledge of the Maxwell and Braginskii stresses in 
order to estimate the rate of th e angular momentum trans- 
port (|Shakura fc SunvaevlHoT^ '). The same problem of the 
magnetic-fleld structure arises in considerations of the ICM 
dynamics and magnetogenesis (see sections 12. II and I7.3|l . 

7.3 Galaxy clusters 

In section [21 we discussed at length the basic properties 
of the galaxy cluster plasmas, the inevitability of pressure 
anisotropics and, therefore, plasma instabilities arising in a 
turbulent ICM, as well as the fundamental theoretical ques- 
tions that this poses. These will not flnd their flnal resolu- 
tion in this paper because it has only analyzed one of sev- 
eral plasma instabilities that must be understood. However, 
just like in the case of the solar wind and the accretion 
flows, an impatient astrophysicist can conceivably glimpse 
the contours of the eventual theory by constraining pressure 
anisotropics and possibly also heat fluxes by the marginal 
stability conditions of the plasma instabilities — with all the 
caveats and uncertainties already discussed above. 

Let us discuss how far this approach can take us in 
answering the three classes of physical problems that were 
described at the beginning of section [^TTI 



7. 3. 1 Regulation of cooling flows 

The apparent refusal of t he galaxy clu ster cores to exhibit a 
cooling catastrophe (e.g.. [Peterson fc"F abian 20])^) has long 
evaded a satisfactory theoretical explanation. A comprehen- 
sive review of the relevant literature is outside the scope 
of this brief discussion. It is probably fair to summarize 
the two main physical mechanisms invoked to explain the 
relatively weak drop in the ICM temperature between the 
bulk and the core as thermal conduction and some form of 
viscous conversion into heat of the mechanical energy in- 
jected into the ICM by the central active galactic nuclei 
(probably in a sclf-rcgulat ing way; se c, e.g., 'Binncv' '20031; 
Kaiser fc Binnov.,2003:,Oinma fc Binncv. 2004 : , Ogrca n ct al.i 
2010l : lTevssier et al.|[2O10r and the references in section [2711) . 
It is clear that the latter mechanism cannot be ignored be- 
cause the thermal condu ctivity of the ICM is unlikely to 
be sufficiently large ('e.g.. FVoigt fc Fabian|[20o3 ) and at any 



rate, thermal conduction is a th ermally unstable m echanism 
of balancing radiative cooling. iKunz et al.l (|201ll ) recently 
proposed that if a sufficient amount of turbulent power is 
assumed to be available, the viscous heating, regulated by 
the pressure anisotropy and, therefore, by the marginal sta- 
bility of the mirror and/or firehose instabilities, can balance 
the cooling in a thermally stable way. They also found that 
assuming such a balance leads to reasonable predictions of 
the magnetic field strength, magnitude of the turbulent ve- 
locities and the outer scale of the turbulence in the ICM. 



7.3.2 Temperature fluctuations and the GTI 

While detailed simulations of the turbulent ICM, bubble dy- 
namics etc. similar to those of lSharma et al.l (|2006l . [ioO?! ) for 
accretion fiows have not been attempted, the marginal sta- 
bility condition for the GTI (equation (|104|l ') could perhaps 
be used to impose a lower bound on the typical scale of tem - 
perature fluctuations in the ICM jSchekochihin et al.ll201(]| '). 
Indeed, if the magnitude of the ion heat flux is limited so 
as to prevent the GTI from being unstable (see section [574} , 
then from equations p04p and poip . we get 



It > P'/^Xmip ~ 5 X 10" 



•,5/2 



1/2 D 

n. B 



1.4 X lO-"' cm, 



(124) 



where Ui is in cm^'^, Ti is in K, B is in G, and the nu- 
merical value has been computed for the plasma parame- 
ters in the core of Hydra A discussed in section 12.11 (see 
equations ((6)1 and (jSj). Interestingly, kpc-scale tempera- 
ture fluctuations are in d eed observed in co o l-core clusters 
llSimioiiescu et al.1 l200l[: [Fabian et al.[ [20061: [Sanders et al.1 



l2010al : [Lagana. Andrade-Santos fc Lima Netd [20101 ). Fur 
thermore, if we use in equation (|124p the physical parameters 
appropriate for the bulk of the cluster plasma, rather than 
the cores (say, Ti ~ 10* K and Ui ~ 10~^ cm""^) we would 
get much larger scales — in the 100 kpc range, which is also 
consi stent with reported obse rvational values for the cluster 
bulk (|Markevitch et al.ll2003l '). 



7. 3. 3 Magnetogenesis 

In the presence of turbulence, the small-scale (fluc- 
tuation) dynamo mechanism generates a mag- 
netic field — this is certainly true in an MHB 
fiuid (e.g.. [ Subramanian. Shukurov fc HaugerJ [200(]| : 
[Brandenburg fc Nordlundl |2009| ')~ How this mechanism 
works in a plasma susceptible to the microscale plasma 
instabilities remains a c ompletely open problem . A rat her 
speculative attempt by [Schekochihin fc CowlevI (|2006h to 
leapfrog the detailed microphysical derivations and model 
the large-scale dynamics based on the idea that the insta- 
bilities would always isotropize pressure towards marginal 
stability values led to a rather dramatic conclusion that the 
ICM might support self-accelerating, explosive dynamos. 
While this conclusion remains to be tested by more rig- 
orous analytical approaches, it does illustrate the general 
conjecture that plasma instabilities are likely to result in 
radical changes of, rather than merely small corrections to, 
the large-scale dynamics of cosmic plasmas. A particular 
mystery in understanding the origin and structure of the 
magnetic field in the ICM is what determines the typical 
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spatial scale of magnetic fluctuations, which observations 
suggest may be substantially smaller than the scale of the 
turbulent motions (see further discussion and references in 
ISchekochihin fc Cowlevii2006i lFl 



8 CONCLUSION 

Let us recapitulate what this paper has achieved and how 
it relates to what was known previously. It has been appre- 
ciated for some time that macroscale turbulence of magne- 
tized weakly coUisional plasma (exemplified by the ICM) will 
naturally produce pressure anisotropics, which will in turn 
trigger firehose an d mirror instabilities at spatial and tem- 
poral microscales l|Hall fc Sciamal 1 19791 : ISchekochihin et alj 
l2005l . see extended discussion in sections [T] and [5}. Since 
the pressure anisotropics are essentially due to local tempo- 
ral change of the magnetic field strength, it is qualitatively 
intuitive that the nonlinear evolution of the instabilities is 
governed by the tendency to cancel this change on aver- 
age; hence it follows that in a driven system (see discussion 
in section [S]) the fluctuations must continue growing in the 
nonlinear regime, albeit s ecularly rather than exponentially 
l|Schekochihin et al.ll2008l . see section \2A\i . 

In this paper, we have constructed a full ab initio 
(weakly) nonlinear kinetic theory of this process for the 
parallel {k± — 0) firehose instability, which is the sim- 
plest analytically tractable case. The evolution not only of 
the fluctuation energy, but also of the full spectrum of the 
resulting firehose turbulence has been worked out, includ- 
ing the effect of gradual spreading of the fiuctuations to 
ever larger scales as the nonlinearly compensated pressure 
anisotropy approaches its marginal-stability value (section 
14. 3[) . We have also extended our kinetic calculation to in- 
clude the effect of ion temperature gradients parallel to the 
magnetic field (parallel heat fluxes). As was pointed out re- 
cently, they lead to a new instability, the GT I, of parallel 
Alfvenic fluctutions (|Schekochihin et al.ll2O10l . see also sec- 
tion [5]2]). Here we have constructed a nonlinear theory of 
its evolution, featuring again a secular growth of magnetic 
fluctuations, but this time developing a spectrum heavily 
dominated by a particular scale (section l5.3p . 

While a speculative discussion of the implications of 
these results for transport in a general magnetized plasma 
(sections I4.4H5.4|) and for particular astrophysical systems 
(section [7]) is possible, a full transport theory has to await, 
at the very least, the completion of similar ab initio kinetic 
investigat ions of the nonlinear evolution of th e mirror in- 
stability (|Rincon. Schekochihin fc CowlevI 120101 1 and of the 
oblique {kj_ 7^ 0) firehosel3 Only then can one attempt 
to devise an effective mean field theory for the macroscale 
dynamics of cosmic plasmas based on solid microphysical 
foundations. A goal of this paper has been to establish a 
template for building these microphysical foundations. 

In the meanwhile, it appears sensible to rely on (or at 



^'^ An example of such observations is IVogt fc Enfilin I |200l), 
although a more recent paper by th e same group appears to revise 
this result llKuchar fc EnfflinllaOOgI) . 

There is a distinct possibility that constructing the most gen- 
eral theory will involve having to study how mirror and fire- 
hose/GTI fluctuations coexist (see section [5]4]l . 



least consider reasonable) the semiquantitative closure ap- 
proach to the macroscale dynamics based on the assump- 
tion that average pressure anisotropies and, probably, also 
heat fiuxes, are set by the marginal stability conditions of 
the microscale plasma instabilities — an approach that has 
found stro ng observational support in the solar wind mea- 
surements (iGarv et al.ll200ll:lKasper. Lazarus fc Garvll2002l: 



Marsch. Ao fc T u"2004V Hellinger et al.ll2006l : lMatteini et all 
2OO7I : I Bale et al.. ,2009. 1 and has already yielded nontrivial 



and possibly sensible physical predictions for the ev olution 
of cosmic mag netism ("Schekochihin fc Co w_le^l2006h . accre- 
tion disk dynamics (Sharma ct al. 2006, 2003), and the tur- 
bulence and heating in the intracluster medium (|Lvutikovl 
l2007l :l Kunz et al.ll201ll l (see further discussion in section [7|. 
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APPENDIX A: KINETIC THEORY: DETAILED DERIVATION 
Al Electrons 

A 1.1 Mass-ratio ordering 

The kinetic equation (|32[) for electrons is (recall that v is the peculiar velocity) 

e /„ Ue X B V X B\ due _ _ 

LE H \ + -T— + Ue ■ V-Ue + V ■ V-Ue 

me \ c c j at 



^ - C[f.], 

9^ (Al) 



where we have labeled all terms according to their ordering in powers of (me/m;)^'^^, while taking Te ~ Ti The ordering has 
been done relative to kvthife and we have assumed 

d T7 ; -1 f meY^^ eB ( mA'^^^ Ue x B 

V"^'/ meCUthe \mij c 

-1/2 

~ OJ, (A3) 



where Vei and i^ee are the electron-ion and electron-electron collision frequencies (they determine the ordering of the collision 
integral on the right-hand side of equation HAip ). We stress that these are formal orderings with respect to the mass-ratio 
expansion, not statements about the exact size of various quantities and their derivatives: thus, some of the quantities ordered 
as unity within the mass-ratio expansion (e.g., kpi or Ue/t'thi) will be ordered small in the subsidiary e expansion to be used 
in solving the ion kinetics (see section 13. 4|) . 

We now expand the electron distribution function in powers of (rrie/mi)^''^: /e — /i"'' + /i^' + • • ■. It turns out that we 
can learn all we need to know from just the two lowest orders in the expansion of equation (|A1|) . Note that we do not expand 
any of the fields — exact E and B are kept. 

A1.2 Order {me/mi)~^ : gyrotropic electrons 
To this order, equation (|A1|) is 

^ '-^"^ ^-nJJ^^O, (A4) 



me 



c dv 9i9 



where Sle = —eB/meC and i9 is the gyroangle variable. Thus, in this order, we have learned that the lowest-order electron 
distribution function is gyrotropic (does not depend on §). 

Al.3 Order (me/mi)"^''^; Maxwellian electrons 
To this order, equation (|A1|I is 

. . V/i°' -^(e+ ^flii^^ • ^ - = (A5) 



me \ c J dv d-d 

Let us multiply this equation by 1 + In /i"' and integrate over the entire phase space. This gives 

dW^Ci/i^'lln/i") =0 (A6) 



becau se the left-hand side of equation (|A5|l is an exact divergence in the phase space. Let us recall that, according to lBoltzmannI 
l|l872h ff-theorem, 

f^J J dVd'w/eln/e = J J dVd^'vCLfelln/e sC 0, (A7) 

where the inequality becomes equality only for a local Maxwellian distribution (the proof for plasmas can be found in, e.g., 
lLongmire|[l963 ). Therefore, equation (|A6|) implies that /i'"' is a local Maxwellian: 

Since v is peculiar velocity, the mean flow Ue has already been accounted for. Note that the perturbation expansion of /e can 
always be constructed in such a way that rie and Te in equation l|A8|) are the exact density and temperature of the electron 
distribution. 
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AI.4 Isothermal electrons 

More can be learned about the electrons without going to higher orders. Let us now substitute the expression 
into equation HA5|) and gyroaverage this equation, (l/27r) J di9, to eliminate the term containing /i^^: 



SI for /, 



(0) 



v«b 



b ■ V?le 



+ 



where Eu = E ■ b. Since equation (jA9[) must hold for all v, it follows from it that 



£11 = -- 



nb ■ Vn, 



(A9) 

(AlO) 
(All) 

The second equation means that electrons (to lowest order) are isothermal alon g the magnetic-field lines, a standard outcome 
of the mass-ratio expansion (|Snvder fc Hammettll200ll:ISchekochihin et 3.1.1120091) . valid up to parallel sc ales ~ Anifp(mi/me) ' 



b ■ VTe = 0. 



(the electron thermal conduction scale; see, e.g., Lithwick fc Goldreich 200ll : Schekochihin et al.ll2009l ). For our fiducial ICM 



parameters, we have Amtp(mi/me) ' ~ 6 x 10 cm, which is larger than the scale I of the motions that have the highest 
rate of strain (see section 1^711 . For turbulent plasmas, this implies globally isothermal electrons (Te = const) because the field 
lines are stochastic. We will adopt this assumption of globally isothermal electrons in all our calculations. 



A 1.5 Generalized Ohm's law 

Let us again go back to equation (|A5|l . multiply it by iJieV and integrate over the velocity space. The result is the electron 
momentum equation to lowest order in the mass-ratio expansion: 

ene (^E + ^ ' j = -Vpe = -TeVn,, (A12) 

where the electron pressure is isotropic because the distribution is Maxwellian, pe = n^Te, and the gradient only affects 
because Te — const fsection lA1.4[) . Note that equation (|A10|) is simply the parallel part of equation (|A12|) . Equation (|A12|) is 
the generalized Ohm's law, equation (|37p . We have thus arrived at the starting point of the derivation in section [32] 



A2 Ions 

A 2.1 Ordering 

The kinetic equation (|32l) for ions is 

^ + u^-v,n + V- v,n + 

at 



Equil 
Pert. 



Ze I „ Ui X B V X B 
— { E + — + 

mi \ c c 



dui 



— Ui ■ ^Ui — V ■ ^Ui 



(A13) 



where we have labeled all terms according to their ordering in powers of e. The ordering has been done relative to fc||Uthi/i 
using the assumptions explained in section [3^ The first row of orderings in equation (|A13|) applies to the equilibrium (lowest- 
order) quantities and their gradients. The second row gives the lowest order in which perturbed quantities appear in each 
term of the kinetic equation. 



A2.2 Expansion of the Lorentz force 

A particular explanation is in order regarding the ordering and the expansion of the Lorentz force. The Lorentz force is given 
in terms of ni and B by equation (|40p . Expanding this equation in e, we have to three lowest orders 

— (e+ ^ ^ \ ^ _ ZnVnii _ ZnV(no,+n20 _ ZnVns, ^ + ... 

nii \ c / mi {noi + nii + n2i) mi {noi + nu) miuoi ^ " So ' (A14) 

where we have used _b| — (see equation (|59[) '). The ordering of the Lorentz force in equation (|A13|) follows from equation 
(|A14|) . Note that, in order to keep precision, we have to keep perturbed densities in the denominators of the first two terms 
on the right-hand side. However, as promised in section [3751 we will see in section [A2. 61 that nu = 0, so the contributions to 
the Lorentz force will start at order and equation (|A14[) will simplify to read 
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Ze 
rrii 



E- 



Ui X B 



ZTe'V {noi + n2i 
miUoi 



I-uaV||^- + ■ 



(A15) 



In section [A2. 91 we will find that n2i = as well. 



A2.3 Order € ^ : gyrotropic equilibrium 

We now proceed to expand the ion kinetic equation (|A13|) . To lowest order, e~^, we get (cf. section lA 1.2 [I 



Zev X Bo dfoi _ ^ dfpi _ ^ 
mi c dv ' i9i9 ' 



(A16) 



where Qi = ZeBo/rriiC. Thus, the ion equilibrium distribution is gyrotropic. We will express the fact that foi is independent 
of the gyroangle i9 by writing foi as a function of two velocity variables, v = and = v ■ bo. In the derivation that follows 
these variables are more convenient than the perhaps more intuitive pair {v±,v^^). Thus, 



fai = foi{t,r,v,v^^). 

Hence follows an identity that will be useful shortly both for foi and other gyrotropic functions: 

'dfor 



dv V \ dv 



+ bo 



dv\\ 



(A17) 



(A18) 



A24 Order e" 

In the next order, equation ([AISP is 



Zev X Bi dfoi 
dv 



rui c 

where we have again used {Ze/mic){v x Bo) ■ d/dv = —Q.id/d'&. Using equation (|A18|I . we get 



dfii _ 1 Ze vx Bi 
di) Sli rrii c 



dfoi 
dvn 



^ (bo X v_ 



Bo 



dfoi 

dva 



Noticing that bo x v± — dv±/d&, we integrate this equation: 

_ Bi fd,hH\ ^ ,^ . 



(A19) 



(A20) 



(A21) 



where gu is an arbitrary function (the gyrotropic part of the first-order perturbed distribution). 

Thus, all we have learned at this order is the gyroangle dependence of fu. This will be a general feature of our expansion: 
since the gyroangle derivative in equation (|A13|) is the lowest-order term, what we learn about each perturbed distribution 
function fu, f2i, fai, ■ ■ ■ , at the lowest order in which it first appears will always be its dependence on i9. 



A2.5 Order e^ : Maxwelhan equilibrium 
At this order, equation ()A13P is 

Ui X B\ dfo 



V ■ Vfu + —iE 

m 



^ df2i ^ Ze vxBj dfu 

dv d'& rrii c dv 



C[fm] 



(A22) 



ii 



where, using equation (|A21|) and other tricks already employed in the two previous sections, we can express the last term on 
the left-hand side of equation (|A22|l as follows 



7i — Q,i [bo X V 



Bo 



Bo 



dv}. 



+ 



dgii 

dv» 



d 



Bo 



d''foi 

dv}: 



+ Va 



Bo 



dgu 

dvn 



.(A23) 



Collisions have made their first appearance at this order and we can now prove that foi is a Maxwellian. The proof is 
similar to the one for electrons in section IA1.3I we multiply equation (|A22|) by 1 -I- In foi and integrate over the entire phase 
space. All terms on the left-hand side vanish because, to the order at which we are computing them, they are all full derivatives 
with respect to the phase-space variables. Thus, 



d'rd^ vC[f o^]lnfo^=0 



foi = 



noi 



2 I 2 



mi 



(A24) 
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A2.6 Order continued: more information about fu 



The fact that foi is a MaxweUian allows us to uncover three important additional pieces of information. First, from equation 
(jA21[) . we learn that fu is gyrotropic: 



fu = gu{t,r,v,v\\). 



(A25) 



Second, we can now prove that nu = 0, as promised in sections 13.51 and IA2.2I Using equations (jA23|l and (|A24|) in equation 
gyroaveraging this equation, (l/27r) J d-d, and substituting for the Lorentz force the lowest-order expression from 
equation (|A14[) . we get 



V\\ I V||/li + — foi 







VllTlH = 0, 



(A26) 



Toi noi J \ Toi 

where the second equation has been obtained by cancelling U|| in the first equation and integrating it over velocities. We 
have used the shorthand Vy = 6o • V, which henceforth will be employed wherever fast parallel variation of the perturbed 
quantities is involved (for slow parallel gradients, we will continue writing 60 ■ V explicitly to emphasize that bo is curved on 
the large scales). Equation ()A26P implies that we may set 

nu = (A27) 

(q.e.d.; see equation (|62p ) and absorb whatever slow- varying density perturbation may arise into noi. After eliminating nu 
from equation (|A26[) . we get 



V||/i, =0, 



(A28) 



so fu has no small-scale spatial variation at all. Note that this confirms equation (|64[l . which was derived from the ion 
momentum equation in the order (i.e., it is the velocity moment of equation HA22|) ') and restricted the fast spatial variation 
of the first-order pressure tensor. Equation (|A15|) is also now confirmed. 



A2.7 Order e^ continued: gyroangle dependence of f2i 



Finally, we go back to equation (|A22|) to determine the gyroangle dependence of f2i. Since fu does not have a small-scale 
part, the first two terms drop out. Using the fact that foi is a MaxweUian and equation (|A23[I . we integrate equation (|A22[I 
with respect to the gyroangle and get 



where 5121 is the gyrotropic part of f2i (so far arbitrary). 



(A29) 



A2.8 Order e^ : role of equilibrium density and temperature gradients 

At this order, equation (|A13|I becomes, upon substitution of the MaxweUian foi and the lowest-order (e'^) expression for the 
Lorentz force from equation (|A15|I 

— - + V ■ V/oi + V||/2i + — — foi 

at Toi noi 



9/3, Ze f vxB2 dfu ^ vxBj 9/2A 2v^^ {V^^uj^-v^ 

iii ^nr -I I ■ — \ • — — I -I 2 JOi = '-'[/liji 



d-d nii 



dv 



dv 



(A30) 



where we have explicitly enforced the assumption that perturbed quantities have no fast perpendicular spatial dependence. 
Analogously to equation (|A23|) . upon using equations (|A25|) and (|A29I) and noticing that v± = —dipo x v)/d-&, we find that 
the last two terms on the left-hand side of equation (|A30|) are a full gyroangle derivative: 



l2 = ni 



d_ 

d-d 



Bi fdfu\ If Bi\ d^fu] , Bi fdg2^\ 2t;||V||<, b^xv^. 



(A31) 



Bo J ydvl j^' Bo \dv^iJ^ vl^^ 

In view of equation (|A31|I and of the gyroangle independence of fu (equation (|A25p ), the gyroaverage of equation (|A30|I is 
dfu 



dt 



+ v« 



J. , ZT, bo- Vtiq. + V||n2, 
bo ■ V/oi + V||32i + — foi 



Toi 



noi 



= C[fu]. 



Since fu has no fast spatial gradients (equation (|A28|l '). averaging equation (|A32[) over small scales gives 
dfu 



at 



t ■ns , ZTe bo ■ Vnoi J. 
bo ■ V/oi + -7^ foi 

-t oi noi 



C[fi. 



(A32) 



(A33) 
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This equation determines fu purely in terms of the equihbrium density and temperature gradients. Since the time variation 
of the equihbrium is slow, it is clear that fu will converge to a steady solution after a few collision times. Then dfu/dt in 
equation (|A33|) can be neglected and the solution obtained by inverting the linearized collision operator. Since we are not 
interested in exact colhsional transport coefficients here, instead of t he full Landau collision op erator, we will use a very simple 
model one — the Lorentz pitch-angle scattering operator (see, e.g.. iHelander fc Sigmaij|2002h . so equation ()A33[) becomes in 
steady state 



^ ^ ZTA bo ■ Vno. _^fv^ 3\ bo- VTo, 



Toi J noi V«fh,- 2 J Toi 



fo^, (A34) 



where ^ — v^^/v and uu is the collision frequency, whose dependence on v is not important here and is suppressed for simplicity. 
The solution of equation (|A34|) that satisfies equation (|A27|) is 



V\\ 

fu = 

Vi-i 



Toi J noi \Whi 2/ Toi 



W|| / 5\ bo ■ VTo, ^ 

foi = [— 77 — 7f fo^■ (A35) 

V<hi 2/ Toi 



The second, simplified, expression above is obtained by noticing that the temperature and density gradients are, in fact, related 
by the equilibrium pressure balance, equation (|66|) . which was obtained in section [331 from the ion momentum equation in 
the order. It is easily recovered by taking the velocity moment of equation HA30|) and averaging out the small scales. Since 
foi is a Maxwellian, the pressure balance takes the form (previewed in footnote [9]) 

1 + ^)^ + ^-0, (A36) 

J-Oi / noi loi 

whence immediately follows the final expression for fu in equation (|A35|I . 

Let us note two useful properties of the solution (|A35[) . First, fu makes no contribution to the pressure tensor: 

Pii J d^vvfu = 0, (A37) 

a result we promised in section [33] (equation ((Mil)- Second, the derivative of fu with respect to ny is isotropic: 

8fu\ _ 1 fv^ 5^So.VT..^^^^^ (^3g^ 



J ^ m V^L 2/ Toi 

which leads to vanishing of one of the terms in equation (|A31|) . 

We will see in Appendix IA2.121 that fu encodes the ion collisional heat fiux (equation (IA35I) is the standard fo rm of the 
appropriate contribution to the perturbed distribution function; see, e.g., equation (D16) of lSchekochihin et al.|[2009l ). We will 
carry the ion-temperature-gradient effect contained in fu through to the end of this calculation because it will interesting 
and instructive to see how contributions from the ion heat fiux arise in the problem. However, this is not the main effect 
we are after and an impatient reader attempting to follow this derivation may find the following simplification useful. If one 
assumes by fiat that VToi = 0, then Vnoi = as well (from equation (|A36I) ) and in all the calculations that follow one may 
set fu = and foi = const, which substantially reduces the amount of algebra. 

A2.9 Order continued: more information about f2i 

Staying at this order, we can learn more about f2i and f-a- Subtracting equation (|A33|) from equation (|A32|l . we get 
^11 fviiga.-f l^^^/o,:") =0 ^ fi + IZl') V||n2, = 0, (A39) 



Tbi noi J \ Toi 

analogously to equation (|A26|I . We have used the fact that, as follows from equation (|A29p . n2i = / d'^vf2i = / A'^vg2i- 
Similarly to the argument in section [A2. 61 this implies that n2i has no fast spatial variation and so we can set 

n2^ = 0. (A40) 

Equation (|A39|) then implies 

V||g2> = 0, (A41) 

i.e., g2i has no small-scale spatial dependence. Since the first term in equation ([A29|) does not contribute to the second-order 
pressure tensor (because the derivative of fu is a function of v only), we have 

P2i = rui j <i'^vvvg2i, (A42) 

and so, in view of equation (|A4ip . P2i has no small-scale dependence. This confirms equation (|67p . derived in section [3.51 
from the ion momentum equation. Note that the tensor P2i will be needed in the large-scale momentum equation (|71[l . It will 
contain the lowest-order pressure anisotropy. 
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A2.10 Order continued: gyroangle dependence of fsi 

Subtracting equation (|A32|) from equation (|A30|) and using equation (|A29|I . we get 



V ,/0i + — JOi 

loi rioi 



Bo 



dfu 
dv« 



(A43) 



where I2 is given by equation (lASlf) (note that the second derivative of fn vanishes there; see equation (|A38|l '). Using again 
the fact that v± — — 9(bo x v±)/d^, we integrate equation HA43|) and get 



dfu 
Bo V dv 

bo X vj_ 
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Bo V^^ii 



2«||V||<, «||V|,B5^ fdfu 

n JOi + 



Bo 



dv« 



+ 



Toi 



foi 



(A44) 



where gsi is the gyrotropic part of /a; (so far arbitrary) and we have used equation (|A36I) to simplify the terms that contain 
equilibrium gradients. 

We will see in section [A2.14l that we do not need to know either g-a or B2 in order to calculate the third-order ion pressure 
tensor Psi and close the ion momentum equation (|68|) for the firehose perturbations. The only remaining quantity we do need 
is g2i — we will now derive the equation for it by going to next order in the e expansion. 



A2.11 Order e'^ 

At this order, equation (|A13[) is, upon substitution of the Maxwellian foi, gyrotropic fu, and equation (|A15[) for the Lorentz 
force, 
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dt 
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Bo 
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mi noi 
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+ bo 
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-^r~ JOi H 2 JOi 



dfu 
dv\\ 



dfu 



V|l<. ) . ( ^ ) = C[f2^]+C[fu,fu\, 



(A45) 



where d/dt = d/dt + uoi ■ V is the convective derivative (with respect to the large-scale flow), /3i = I'thi/'i^A- Note that V/ii in 
the above equation is with respect to slow spatial variation. Note also that the collision operator at this order has two parts: 
the linearized operator describing interaction of f2i with the Maxwellian equilibrium foi and the nonlinear operator, denoted 
C[fu,fu], describing interaction of fu with itself. 

We will only ever need the gyroaverage of equation (|A45p . Many terms then vanish or simplify. What happens is mostly 
straightforward: the gyrovaerages of d/d-ff are zero, the gyroaverages of the velocities are done using the identities 



(v) = U||bo, (vv) = — ( ' ^ t>obo 



- viibobo 



(A46) 



(henceforth angle brackets denote (l/27r) J di9). There are a few terms that are perhaps not obvious and so require explanation. 

First consider the term 1; ■ V/ii. We showed above that fu is a function of v, «||, and r. However, the spatial gradient 
here is still taken at constant v. Since one of the new velocity variables = d • bo is a function of v and r, we have 



« • (V/iO„ = V ■ iVfu),^,,^ ^[vv: Vbo 



dfu 
dvu 



('^ • (V,fiO J = «||bo ■ V/n + ( V ■ So 



dfu 

dv« 



(A47) 



In the final expression, V/i; is now understood to be at constant v and Since V ■ So = — ^So • Vi?oj /Bo, the additional 
term that has emerged is readily interpreted as the mirror force associated with the large-scale variation of the magnetic field. 
Now let us turn to the two terms in equation (|A45|) denoted by I3: the second of these terms gives 



Zev X B2 df2, 
rrii c dv 
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Bo 
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+ bo 
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(A48) 



where we have used equation (|A29|I and the fact that (9/ii/9u||)^ only depends on v (equation HA38fl '): the first term, upon 
substitution of equation (|A44p . gives 



' Ze vx Bj dfM 
rrii c dv 



Zev X Bi 
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Bo 
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Bo 
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(A49) 



The first term in equation (|A49|) exactly cancels when equation HA48|I is added to it. Simplifying the double vector product 
in the second term and noticing that the gyroaverage in the third term is equal to ^1 — 5oSo^ , we have 
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2 At Bl + 



Bo 



4 Bl 



dfu 

dv\\ 



+ v 



+ v 




(A50) 



where we have used the perturbed induction equation (|60p to express Vyitii in terms of Bi . Finally, the gyroaverage of 
equation (|A45|I is 
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(A51) 



where (/s) is given by equation (|A50[) . Note that g2i does not have fast spatial variation (equation (|A41|) '). so, to the order 
we are keeping, dg2i/dt — dg2i/dt. The next step is to average this equation over small scales: again many terms vanish (in 
particular, all terms where the fast-varying perturbed quantities enter linearly) and we get 
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/o> = C[52.] + C[/ii,/ul, 



(A52) 



where the overline denotes the small-scale average and the derivatives of fu have been written in a compact form that will 
prove useful momentarily. Remarkably, the contribution of the perturbations has survived in equation (|A52[) in the form of 
a single quadratic term — in section IA2.13I we will see that it gives rise to precisely the nonlinear feedback on the pressure 
anisotropy that was anticipated qualitatively in section 12.41 

By averaging out the gyroangle- and small-scale-dependent parts of the third-order kinetic equation HA45|) . we have 
eliminated fa, <73i, B2, -B3, and U2i, which are unknown and potentially very cumbersome to calculate. As we are about 
to see, in order to calculate the ion pressure tensor to the relevant orders, we do not, in fact, need to know any of these 
quantities, so we will neither have to revisit the unaveraged equations (|A45P or (|A5H) or go to higher orders in the e expansion 
of equation (|A13|I . Equation (|A52|I determines g2i, which is all that we require to calculate P2i (section lA2.13|l and Pa; (section 
IA2.14|I . Knowing these tensors will then allow us to close the ion momentum equations describing the plasma motion at large 
(equation ([71}) and small (equation (|68j) scales. 



A2.12 Order e'^ continued: transport equations 
Since 5121 is gyrotropic, the tensor P2i is diagonal: 



P2, = P2i ( I - bobo ) + P2,bobo = P2i\ + (P2i - P2J ( 3 ' ^ ''"''0 



(A53) 
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where the scalar pressures are 

± /■ ,3 ruivl II /■ 3 2 2 X 1 II f ,3 rriiv'^ 

P2i = d V g2i, p'^^= dvmiV^\g2i, V2i ^ ^P2i + ^vli = <^ v g2i. (A54) 



Then the ion momentum equation (|7H) become; 
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duoi „ - _ 
minoi—— = - Vp+ V 
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bobo ( P^, - pi; + ^ 



P = P2. + i(P2^,-pL) + ^, (A55) 



which is the famihar momentum equation in the long- wavelength limit (equation H18|) '). Let us first explain how p2i (or, 
equivalently, p) is determined and then calculate the pressure anisotropy — p^^ (section |A2.13|) . 

First, let us integrate equation (|A52|) over velocities. Since J d'^v g2i = and J d'^vv^^fu = 0, we get 

= -noiV ■ uoi, (A56) 

dt 

an unsurprising result (the continuity equation was already obtained in section [331 see equation (|70|) ). Now multiply equation 
HA52|) by niiV^ /3 and integrate over velocities: 

dpoi , ^P2^ 2 /- \ 5 

^ + ^ = "3 ^ ■ r"''') - 3"°'^ ■ 

where poi = noiToi and the parallel coUisional heat flux is, using equation (IMsliFI 
3 m^v'^vu /■ 3 m,v^vl f 5\ So ■ VTo, 5 



/ d%—-^fu = - <iv—-^l—^^]fo. " ' ^-^no,^bo- VTo,. (A58) 

The heat flux term in equation (|A57I) arises from the third and fourth terms on the left-hand side of equation ()A52p (in the 
fourth term, write v\ — — and integrate by parts with respect to at constant v). 

Since formally all terms in equation HA57P except the one involving p2i have slow time dependence, we may assume 
that so does p2i and, therefore, dp2i/dt can be dropped from this equation (this can be formalized via averaging over short 
timescales). Using now equations HA56|) and (|A58|) . we can rewrite equation (|A57|) as an evolution equation for the equilibrium 
temperature: 

noi^^^ = V • fnoi^iSoSo • VToi) — ^ noiToiV ■ uo,, Ki = ^-^^i (A59) 
dt \ J 6 6 Vii 

where is the ion thermal conductivity. The first term on the right-hand side represents coUisional heat transport, the second 
compressional heating. 

Equations ()A56|I and (|A59|I evolve noi and To;. However, the equilibrium density and temperature can only change in 
such a way that pressure balance, equation HA36|) . is maintained. This means that if we know the spatial distribution of Toi, 
we also know that of noi , or vice versa. Compressive motions will develop to make the density and temperature distributions 
adjust to each other and preserve the pressure balance. These motions must be consistent with the momentum equation 
l|A55|l and the isotropic pressure perturbation p2i will adjust to make it so. Thus, if we provide the expression for the pressure 
anisotropy P2i — p\i, the other 5 equilibrium quantities — itoi, noi, Toi, P2i, and Bo — are determined by the closed set 
of 5 equations momentum equation HA55|) . continuity equation HA56p . heat conduction equation (IA59|I . pressure balance 
equation (|A36|I . and induction equation (|69p . 



A2.13 Order continued: pressure anisotropy 

In order to calculate the pressure anisotropy P2i — p|i, we multiply equation HA52|I by niiiv^ — 2uj|)/2 and integrate over 
velocities. All isotropic terms vanish, as does the term containing dfn/dv, which generally cannot contribute to pressure 
(after integration by parts, it is zero because J d'^vvfn = 0). The result is 

^ (P2i -pL) = 3po, f SoSo : Vuoi - i V • UQi 4- ^^ ^"^^J j - V ■ [So(gu - <?l'i)] - 3qnV ■ bo - Sm(p2i -pL), (A60) 



Formally, this equation is the result of taking the velocity moment of the kinetic equation IIA13II at the order t^. We do not write 
explicitly equation l)A13|l at this order because it is not needed for anything except the momentum equation, the form of which we 
already know. 

•^^ Since we used a very simplified collision operator in our calculation of fu in section IA2.8I the numerical prefactor in the expression 
for the heat flux should not be regarded as quantitativ ely correct. This i s not a pr oblem for our purposes. Correct numerical coefficients 
for this and other coUisional fluxes were calculated by iBraginskiil l l 19651) (see also lCatto fc Simakovll2004l) . The same caveat applies to 
equation IIA61II and equation I IA65I I. 

If, as discussed in section [A2.8I one takes the easy option and assumes VToi = and Vnoi = 0, then no equations are needed for 
the constant density and temperature, while p2i in equation l|A55[l is determined by the incompressibility condition V ■ uoi = (which 
follows from equation I IA56II with noi = const). 
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where q^^ and ql^ are parallel fluxes of perpendicular and parallel heat, respectively: using equation (|A35|I 

gii = / d w fu = --no, 6o • vToi, q"^ = d v m^v^^ fu,= "-Oi oo ■ vTo,, (A61) 

and we note that q^^ + Qii/Z = gii (see equation ([X58])). To work out the collision term in equation (|A60|I . we have again 
resorted to brutal simplification by using the Lorentz operator (see equation HA34P ') and dropping the nonlinear collision term 
C[/ii, fii] in equation ( A52I) (the consequences of retaining this term, which are mostly small and irrelevant for our purposes, 
Catto fc Simakovlliooi 'l. The solution of equation (|A60|) is 



have been explored by 



A(t)^^^Hi^ = Ao + | Tdt'e— (A62) 

POi Jo of t)Q 

where Ao is the part of the anisotropy due to the large-scale dynamics: 

if., 1 V ■ \bo{qh - q\S\ + 3gH V ■ b„ ] 

Ao = — <^ bobo : Vno, - - V ■ uo, 5-^ > (A63) 

Vii 3 3poi I 

and we have assumed A(0) = Ap. The first two ter ms in equation (|A63|I are the well known coUisional contributions to the 
pressure anisotropy calculated by iBraginskiil (j 19651 ). The heat-fiux terms did not occur in Braginskii's calculation because 
they were small in his assumed sonic-fiow ordering (uoi ~ t'thi)- They occur here because our ordering is subsonic (uoi ~ 
efthi; see equation (I52t) — that heat fluxes appear in the press ure tensor under such assumptions is also a known fact 
l|Mikhailovskii fc Tsvpinlll97ll . 1 1984 ICatto fc Simakovll2004l . I2OO5I I . The new, nonlinear part of the anisotropy is the second 
term in equation HA62|I . which is due to the flrehose fluctuations. This result was predicted on heuristic grounds in section 

[m 

Equation (|A62|) completes the set of transport equations derived in section IA2.121 but we still need to calculate |Bj^p. 
This is done via equations (|60l) and (|68[) . The third-order pressure term in the latter equation is be calculated in the next 
section. 



A2.I4 Pressure tensor for the firehose turbulence 



In order to close the small-scale momentum equation (|68p . we must calculate the divergence of the third-order ion-pressure 
tensor or, more precisely, the perpendicular part thereof: since the fast spatial variation is only in the parallel direction, we 
have, from equation (|A44[) 



dg2i 
dv« 



Bi 2 bo 



~?2 JOi + 



Bo 



dv\\ 



poiVyiiii - {2qu - qii 



V||gi 

Bo 



X 60 >, 



(A64) 



where the last formula was obtained via integration by parts with respect to uy (at constant i,). Note that the terms in 
equation (|A44|) containing B2 and gzi do not contribute, so, as announced at the end of section IA2.11I we do not need to 
compute these quantities. The pressure anisotropy — is given by equation (|A62|) and the heat fluxes q^i and by 
equation (|A61I) . whencj^ 

^ _ 2qi, - gll, 1 TJth, So ■ VTo, 



- o rr. ■ (A65) 

POi^'thi 2 Vu Toi 

The second term in equation (|A64|I is recognizable as the collision-independent "gyro viscosity" j Braginskii' '1965') a nd the 
third term as the coUis ional heat-flux contribution to it that arises for subsonic flows (|Mikhailovskii fc Tsvpin 1971i . Il984l : 
ICatto fc Simakovll2004h . It is the gyroviscous term that will limit the range of wavenumbers susceptible to the firehose 
instability (see section W7I\ . while the heat-flux term will lead to substantial modifications of the firehose turbulence and even 
give rise to an additional source of unstable behaviour (the gyrothermal instability; see section [5)|. 



A2.15 CGL equations with nonlinear feedback 

It is perhaps useful to explain how our equations compare to the standard ones, due to lChew. Goldberger fc Low! (|l956l ). or 
CGL. Let us first notice that the induction equation (|69|l implies 



2* As explained in footnote 1261 the numerical prefacto r here should not be take n literally because we have used a very simplified collision 
operator. The correct prefactors can be found, e.g., in lCatto fc Simakovl | |2004|) . 



34 M. S. Rosin, A. A. Schekochihin, F. Rincon and S. C. Cowley 



-^-rr ^ ^obu : Vnoi - V • uoi. (A66) 
Using equation HA56|I for V ■ uoi, we may, therefore, rewrite equation HA60|) as follows 
d , A_ II , ^ f 1 dB 2 1 dnoi 



where B includes both the large-scale magnetic field Bq and the averaged firehose fluctuations (see equation (|28p ). Now recall 
that the total perpendicular and parallel pressures may be written as follows: 

1 2 I 

= POi +P21 + - (pi -pli), p1' = POi +P2i- - (pii -pli), (A68) 

and so, combining equations (|A57|I and HA67|) . we get 



— In = = ^ V,, (A69) 

at noiB Poi Poi 



± In ^ = ^ 1 2m ^' ^' . (A70) 

at ngj poi Poi 

These are exactly the CGL equations (without n eglecting the heat fluxes; see also ISnvder. Hammett fc Dorlandl Il997l : 
ISnvder fc HammettI [iooH : iRamodlioOsI : iPassot fc Sulcm, 20071 . who adapt these equations to a situation in which collisions 
are weak and the heat fluxes a re calculated from wave-pa rticle interactions). They are also the equations that form the basis 
for recent numerical studies b ylSharma et all (|2006l, |2007| ) of the effects of pressure anisotropics in astrophysical plasmas. As 
pointed out in ISharma et al.l (|2006l ). equation (|A69[) can be rewritten in a form that makes explicit the conservation of the 
first adiabatic invariant (including by the heat-flux terms): 

d pi (pi - qi \ pi -p! , . s 

+ V ■ =- uoi + bo =- = -Vii -. (A71) 

dt B \B B J Poi 

T he non linear feedback in equations (jA69f[A70ll is provided by the small-scale magnetic fluctuations in _B. ISharma et al.] 

l|2006l . |2007| ) do not have the small-scale fluctuations and model their effect by introducing strong effective damping terms 
in equations (|A69flX70)) that limit the pressure anisotropics to the marginal state of the plasma instabilities. The calculation 
carried out in the present paper attempts to derive this feedback from first principles (note the discussion in section 16.11 
regarding the absence in our theory of particle scattering by firehose fluctuations). 

A2.16 Comparison with the equations of \Shaviro & ShevchenkJi 1 196 A) 

Finally, l et us make a comparison between our equations and those derived in the classic paper by IShapiro fc Shevchenkd 



l| 19641 ') n Their approach is to assume isotropic electrons, an initial bi-Maxwellian equilibrium distribution for the ions with 
p±_ 7^ P|| , no collisions, a constant uniform background magnetic field Bq, and to calculate the evolution of the ion equilibrium 
in the quasilinear approximation. Recast in our notation, their equations (12-12') for the parallel (fc^ = 0) firehose read 



dp,^ _ ± d 1 15BxP „ / X ||\ d 1 |<5BxP d \5B 



At dt2 Bl "V" '^Vdt2 Bl At Stt ' '^'^''^^ 



At - ^^^At2 Bl ^^)At2 Bl ' 

The first terms on the right-hand side of these equations are the CGL terms and by t ransferring them to the left-ha nd side, 
we recover equations (jA69|l and (|A70[) with zeros on the right-hand side. Indeed, since IShapiro fc Shevchenkd (|l964 ) assume 
Vii = 0, Bo = const, and bi-Maxwellian ions, there is no collisional relaxation, dln_B/dt = (l/2_Bo)d|5BxP/di (see equation 
(|28p '). and there are no flows or heat fluxes. The second terms in equations ()A72[I and (|A73|I could be interpreted as the 
relaxation of the pressure anisotropy by effective scattering of particles off the firehose fluctuations. However, these terms, as 
well as the last term in equation (|A72|) . are subdominant in our ordering (~ compared to the C GL terms; see equation 
(|53p and the discussion in section [STT]) . They are, in fact, also negligible under the assumptions that IShapiro fc Shevchenkd 
(|1964 ) have to make in order to guarantee the validity of the quasilinear approximation, viz., pf — pi <C p| and /3i >• 1 (see 
discussion at the end of their paper). 

Note that, dropping the subdominant terms and so retaining only the CGL terms, one can integrate either equation (|A72[) 
or equation (|A73P and get equation (|12ip . which then leads to equation p20p if A(t) is set to its ma r ginal value. This is the 
result for the saturated fluctuation amplitude expressed by equation (17) of IShapiro fc Shevchenkd ( 19641 ). As discussed in 



section [6^21 such a constant saturation level is only possible in a system that is coUisionless and where anisotropy is assumed 
to relax from an initial value rather than continuously driven. 

We thank the anonymous referee for suggesting to include this discussion. 



